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(3JT)' We derive a set of coupled four-dimensional integral equations for the 

NN — ttNN system using our modified version of the Taylor method of 
classification-of-diagrams. These equations are covariant, obey two and three- 
body unitarity and contain subtraction terms which eliminate the double- 
counting present in some previous four-dimensional NN — ttNN equations. 
The equations are then recast into a from convenient for computation by 
■ grouping the subtraction terms together and obtaining a set of two-fragment 

scattering equations for the amplitudes of interest. 

O 
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^ ! I. INTRODUCTION 

OS f 

The iViV — TtNN problem occupies a privileged place in nuclear physics. Not only does 
its history stretch back to Yukawa's original efforts to model the strong nuclear force [|T|, 
but it is still the subject of considerable research to this day. This current interest occurs 
partly because the iViV — ttNN system is one light nuclear system which can be used to 
test the validity of models in which hadrons are the fundamental degrees of freedom. Only 
in systems which, like the iViV — ttNN system, contain a relatively small number of degrees 
of freedom, is it possible to complete a calculation of experimentally observed quantities 
within the framework of such a model, while retaining some control over the approximations 
used. Consequently, considerable effort has been put into attempting to produce physically 
accurate hadronic models of the NN — ttNN system. The predictions of these models can 
then be compared with the mass of available experimental data and a judgement on the 
models' validity formed. 

Over the past twenty years considerable theoretical progress has been made in this di- 
rection, with the culmination being the independent derivation, by a number of groups 
using different techniques, of a set of scattering equations known as the NN — ttNN equa- 
tions 0d3- ^he P r °bl em is that models based on these NN — nNN equations which 
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treat the re — N amplitude as the sum of a pole and non-pole term fail to reproduce the 
experimental data. This may well be because of certain theoretical inconsistencies in the 



equations. For example, Jennings [|13|] has pointed out that the NN — tcNN equations in- 
clude the diagram on the right of Figure [I], but exclude that on the left — even though the 
left-hand diagram is merely a different time-order of the right-hand one and both diagrams 



represent the same physical process. Jennings and Rinat [13,14| and Mizutani et al. [15 
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have made it plausible that the omission of this diagram accounts for the failure of models 
based on the NN — ttNN equations to correctly predict the tensor polarization, T 20 for -rr — d 
scattering. (For more detail on the history of the theory of the NN — ttNN system and 
a thorough comparison with the experimental data, see the recent book by Garcilazo and 
Mizutani ||I§.) 

The standard NN —ttNN equations do not include the Jennings mechanism because they 
are derived by using unitarity as a criterion for truncating the full field theory of nucleons 
and pions. One way in which this derivation can be done is to examine the diagrammatic 
expansion for the NN — ttNN amplitudes in old-fashioned or time-ordered perturbation 
theory (TOPT). The NN — nNN equations may then be obtained by truncating this ex- 
pansion at the one explicit pion level. (See [11],|17[] for details.) This truncation not only 



appears to lead to incorrect predictions for T 2 o, but also produces the under-dressing of 
the two-nucleon propagator discussed in |TJ,[n|. Indeed, the fundamental shortcoming of 
the standard iViV — nNN equations is that two diagrams each representing the same set 
of physical processes, but having them occur in a different time-order may end up treated 
completely differently. In fact, as we see in Fig. ^ one may be included and one excluded, 
with potentially disastrous consequences. 

This difficulty can be entirely circumvented, however, if the diagrammatic expansion of 
the field theory of nucleons and pions is performed in covariant perturbation theory, rather 
than time-ordered perturbation theory. This is so because TOPT diagrams which differ only 
in the time-order of the physical interactions involved are all included in the one covariant 
perturbation theory graph. In particular, the two TOPT diagrams in Figure |l] arise from the 
same covariant perturbation theory graph. Hence, if equations for the NN — nNN system 
are derived from this covariant perturbation theory diagrammatic expansion they will not 
suffer from the fundamental deficiency of the standard iViV — nNN equations. 

Pursuing this approach to the NN — nNN system also has the advantage that, provided 
care is taken in truncating the field theory, the equations will automatically be covariant, 
since they are derived from a covariant diagrammatic expansion. Furthermore, if inte- 
gral equations are used to sum classes of diagrams containing infinitely many perturbation 
graphs, they generate amplitudes which are non-perturbative, even though the equations 
themselves may originally be derived from the perturbative diagrammatic expansion. Ex- 
amples of this are the Schwinger-Dyson equations of a field theory which may be derived from 
a resummation of the original perturbation expansion of the theory, but in no way involve 
a truncation of this expansion at some order in the coupling constant of the theory. Such 
integral equations governing the NN — nNN system will of necessity be four- dimensional, 
and hence the numerical solution of these equations in order to obtain experimental predic- 
tions is a challenging problem. However, until this problem is tackled it is not clear that the 
one-explicit-pion sector of the NN — tiNN system has been properly dealt with. Indeed, 
the standard theory of the iViV — wNN system, based on time-ordered perturbation theory, 
only approximates the fuller description based on a covariant diagrammatic expansion. Un- 
til it is clear that this approximation is an adequate one appeals to mechanisms beyond the 
one-pion sector to remedy the present disagreement between theory and experiment cannot 
be definitely upheld or overturned. The derivation of equations for the iViV — nNN system 
in the framework of covariant perturbation theory is therefore an important question at the 
heart of one of the oldest problems in theoretical nuclear physics. 
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The question then is: how are such equations to be derived? Given a Lagrangian one 
could use functional techniques to derive the relevant Schwinger-Dyson equations of the field 
theory, truncating the Schwinger-Dyson equation hierarchy by some approximation scheme, 
and so obtaining a set of four-dimensional coupled integral equations for the NN — nNN — 
-kitNN — . . . system. In this paper we choose not to employ such a functional calculus 
technique, but instead use the Taylor method of classification-of-diagrams |20|. 

The Taylor method is a general one allowing the derivation of equations connecting 
the amplitudes obtained from diagrammatic expansions, and hence is ideally suited to the 
derivation of equations from the Feynman diagrammatic expansion of the amplitudes for 
the iViV — nNN system. Reference |2T|, henceforward known as paper I, presented a review 
of the Taylor method. That paper also pointed out that Taylor's method leads to double- 
counting when applied to certain covariant perturbation theory amplitudes. It was shown 
how the Taylor method could be modified in order to eliminate this double-counting, thus 
producing a technique by which double-counting-free covariant four-dimensional scattering 
equations may be derived for, not only the NN — nNN system, but also other few-hadron 
systems, such as the nN — irnN system and the problem of pion photoproduction. 

This idea of applying Taylor's original method to the NN—tcNN system is not a new one. 
Three pairs of authors have already attempted to derive covariant four-dimensional equations 
for the NN — nNN system using Taylor's method. Firstly, Avishai and Mizutani (AM), 
derived coupled covariant four-dimensional equations for the iViV — nNN system using the 
Taylor method |IJ . However, Avishai and Mizutani failed to eliminate the double-counting 



which arises when the Taylor method is applied directly to the graphical expansion of a time- 
dependent perturbation theory. Consequently the equations they derived double-counted 
certain diagrams. Avishai and Mizutani then used Blankenbecler-Sugar [22] reduction in 
order to reduce their covariant four- dimensional integral equations to more manageable 



three-dimensional ones [p3|l . The reduced equations thus obtained are equivalent to the 



standard three-dimensional iViV — itNN equations AM (and others) had previously derived. 
The standard NN — irNN equations do not contain any double-counting and thus the final 
set of equations used for numerical work by AM may be regarded as double-counting free. 
However, AM's "derivation" of these equations by a three-dimensional reduction of four- 
dimensional equations which themselves do contain double-counting is open to question. 
Secondly, in 1985 Afnan and Blankleider (AB) derived a set of covariant BB—nBB equa- 



tions, in which the baryon B could be either a nucleon or a delta [17|. However, instead of 



using Taylor's original classification-of-diagrams scheme Afnan and Blankleider used a sim- 
plification of Taylor's method which they, together with Thomas and Rinat, had developed 
some years before P,|TT| . (This simplification and its relation to Taylor's original work were 
discussed in paper I.) The equations thus obtained by Afnan and Blankleider were exactly 
those found by Avishai and Mizutani, except that the nucleon iV in Avishai and Mizutani's 
theory was replaced throughout by a baryon B which could be either a nucleon or a delta 
[]. The use of a diagrammatic technique similar to that used by Avishai and Mizutani, and 



x There were terms present in AM's equations which were not included in AB's results, but the 
addition of these terms to AB's equations could have been effected with only minor changes to 
their argument. 
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the consequent derivation of similar equations naturally meant that Afnan and Blanklei- 
der's equations also contained double- counting if viewed in a four- dimensional framework. 
However, as in Avishai and Mizutani's case, this problem was never fully revealed, since the 
BB — ttBB equations were solved in a three-dimensional time-ordered perturbation theory, 
and so the double-counting was temporarily hidden, even though it was still present in the 
full four-dimensional theory. 

Recently, Kvinikhidze and Blankleider (KB) have recognized the double-counting in these 



two derivations |2Jj2^|. They have introduced a modification to Taylor's method as applied 
by AM and AB which allows them to derive a set of covariant four-dimensional equations 
for the NN — nNN system which are free from double-counting. These equations are, 
apart from a minor point, equivalent to the ones derived here. However, this work differs 
from that of Ref. [0,0 in three main ways. Firstly, the method used here is rigorously 



based on the modification of Taylor's original classification scheme for anm->n amplitude 
in any perturbation theory as detailed in I. By contrast, KB have used a classification-of- 
diagrams scheme which is only loosely based on Taylor's original work. Secondly, our use 
of the full Taylor method allows us to exploit the true Lagrangian independence of that 
technique. Hence here we do not specify the Lagrangian to be used in the description of the 
NN — nNN system, while KB restricted themselves to the case of a 3 field-theory. Finally, 
we go beyond KB's work in casting our equations in a form convenient for computation by 
deriving a set of coupled equations for two-fragment amplitudes in the NN — nNN system. 

This work proceeds as follows. In Section [H| Green's functions and amplitudes of the 
field theory are defined while Section |HI| provides a brief summary of the Taylor method of 
classification of diagrams. In Section [TV] the iViV amplitude T^n is discussed and is found 
to depend on the fully-connected two-particle irreducible NN — > ttNN amplitude, . 
Section [V] is therefore devoted to deriving an equation for a process which involves 

eliminating the double- counting that occurs in the equation found when the Taylor method 
is applied directly to F^ 2 ^. Section VI then outlines the removal of double-counting from 



the equation for the two-particle irreducible amplitude T^ N . Once in possession of double- 



counting-free equations for the amplitudes involved we proceed in Section [VI I| to derive a 



set of coupled equations for the NN — ttNN system. Section |V111| then explains how these 
equations may be anti-symmetrized in order to obtain the physical amplitudes. Finally, 
Section |IX| gives suggestions for the specification of the input amplitudes and explains how 
physics beyond the one-pion sector can be implicitly included in the framework developed 
in Sections 



The work concludes with four appendices containing details of arguments to do with the 
dressing of the propagators, the amplitudes which are the input to the equations, and details 
of the double-counting removal performed in Section [VT] . 



2 A summary of these results was presented at the 14 th International Conference on Few Body 



Problems in Physics [26|. 
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II. THE FIELD THEORY OF NUCLEONS AND PIONS 



In this section we define the field theory of nucleons and pions used in this work by 
discussing the Lagrangian to be used and defining the Green's functions and amplitudes of 
the theory. 

A. The Lagrangian 

The Taylor method is essentially a topological classification of the diagrams in a per- 
turbation expansion. Therefore, it depends only on the topology of the diagrams which 
can be generated, and not on the details of the particular Lagrangian under consideration. 
Consequently, in this work we do not specify the precise form of the Lagrangian to be used. 
We merely assume that the Lagrangian density describes a system of nucleons and pions 
and so takes the form: 

£ = C'T) + £</> + £-inti (1) 

where: 

C D (x)=^(x)(i-fd li -m)^(x), (2) 
= \(dj(x) ■ P$(x) - m 2 J(x) ■ $(x)), (3) 

— * 

where ip is the nucleon field and 4> the pion field, which is an isovector. 

At this stage we do not actually need to assume anything about the structure of Ci nt 
in order to apply the Taylor method. However, for concreteness the reader may wish to 
assume that Ci n t is a sum of a bare NNn vertex, proportional to ifxpip and a contact term, 
proportional to i/jcfxpip. As will be seen below, our derivation never demands that we say 
more than this. Bare form factors for these two kinds of vertices could be included in the 
Lagrangian if we wished. Such form factors would be a natural place to include information 
from QCD about the structure of the nucleons and pions in our field theory. However, for 
the present we leave these details to one side. 

Note also that the Lagrangian independence of the Taylor method makes it relatively 
easy to return at a later stage and add other mesons and baryons to the field theory. The 
addition of extra particles to the theory merely complicates the resulting equations without 
changing the spirit of the derivation. However, as a first step in this paper we have derived 
the results for a field theory in which the only quanta are nucleons and pions. 

B. Green's functions and amplitudes 

The coordinate space Green's function for the transition from a j'-nucleon, m — j-pion 
state to a j'-nucleon, n — j'-pion state is defined by: 

a( i',i) i i i i i \ 
\J / I rp rp np rp • rp rp . rp . rp \ 

n<— mv^l ' • • • 5 ^j' i ^j'+l i ' ' ' i ^ni ^l? • • • i ^3 1 ^j+l i • • • i ^m) 

((TOK) . . • ■ -tKMxi) ■ ■ .^(xMxj+i) ■ ■ -#O)|0>, (4) 
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where |0) is the vacuum state, T is the time-ordering operator, is the nucleon an- 

nihilation (creation) operator, and <p is the pion operator, whose isospin index has been 
suppressed. 

This Green's function may then be Fourier transformed in order to obtain the momentum- 
space Green's function: 

^n<— miPl i • • • i Pj' i Pj'+li • • • i Pm Pit ' ' • > Pji Pj+1 > ' ' ' ) Pm) ■ (5) 

From this momentum-space Green's function the amplitude: 

■^n*^m{p'li ■ ■ ■ i P'j>iP'j'+ll ■ • • > Pn'i Pi) ■ ■ ■ ! Pj)Pj+li ■ ■ ■ j Pm) (6) 



is obtained via the LSZ reduction 27 formula: 



A&%(j/ lt ■ ■ .,p' n ;pi, ■ ■ ■ ,p m ) = Z N ~ dr^{$>\) ■ ■■dN 1 (p' f )d~ 1 (p' r+1 ) . . . d~VJ 

G^(p'i, ■ ■ ■ • • • ,Pm)d- N \pi) . . .d N \ Pj )d~\p 3+l ) . ..d-\ Pm )z N "z^ (7) 



where and d^ are the free single-nucleon and single-pion propagator in momentum space, 
given by the formulae: 

dN ^ = nr 1 — ' ( 8 ) 

p — 771 TV 

7T 

Note that we are using fully dressed propagators in order to do the reduction here and so 
Z/v and Z n are the wave function renormalizations for the nucleon and pion respectively, 
while rajv and are the dressed nucleon and pion masses. The Green's function G^j^, 
and therefore the amplitude A^J*^, may be expressed as the sum of a perturbation series of 
Feynman diagrams. 

Note that although all the formulae in this section are written for identical particles, 
in this paper we (at first) consider only the amplitudes for distinguishable particles. The 
physically correct indistinguishable-particle amplitudes may then be obtained from these 
distinguishable-particle amplitudes by the usual symmetrization and anti-symmetrization 
processes. This procedure is fully justified and implemented in Section |V111. 



III. THE TAYLOR METHOD OF CLASSIFICATION-OF-DIAGRAMS 

A. Taylor's original method 

Taylor's method of classification-of-diagrams then provides a scheme for classifying all 
of the diagrams which contribute to any amplitude according to their topology. In 

this paper we concentrate on the topological structure of the diagrams in the s-channel 
0, although it is possible to apply the Taylor procedure in other channels instead of, or 



5 We use the notation s—,t— and n-channel throughout this paper, intending it in the sense defined 



by Mandelstam f2|,|2 
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as well as, applying it in this channel. Note that because of this focus on the s-channel 
structure, unless otherwise specified, the irreducibilities given for amplitudes and diagrams 
are s-channel irreducibilities. Note also that from now on we suppress the indication of the 
number of nucleons and pions present in the initial and final states, in order to simplify the 
notation. 

Having noted these points the Taylor scheme then requires the following definitions: 

Definition 1 (r-cut) An r-cut is an arc which separates initial from final states and inter- 
sects exactly r -lines, at least one of which must be an internal line. If all of the r lines cut 
are internal lines then the cut is called an internal r-cut. 

Note that in writing this definition we assume that all perturbation diagrams are rep- 
resented in a two-dimensional plane. We do allow the lines in any diagram to "jump over" 
one another: two lines do not have to meet at an interaction vertex whenever they intersect. 
By contrast, a cut is defined to intersect all the lines it meets: it may not jump over any of 
them. (Other definitions of an r-cut, which do not assume that the diagrams are represented 
in the plane, may be composed but it is the above definition which Taylor himself used.) 

Definition 2 (r-particle irreducibility) A diagram is called r-particle irreducible if, for 
all integers < k < r, no k-cut may be made on it. An amplitude is called r-particle 
irreducible if all diagrams contributing to it are r-particle irreducible. 

The r-cuts which may be made on a particular diagram are now divided as follows: if an 
r-cut is not internal it is called "initial" if it intersects at least one initial-state but no 
final-state line; "final" if it intersects at least one final-state, but no initial-state line, and 
"mixed" if it intersects both initial and final-state lines. 

Once these definitions are made any diagram contributing to the connected s-channel 
(r — l)-particle irreducible (r — 1PI) m — > n amplitude, Affi^l, may be placed in one of the 
five classes C\ to C5. The class a particular graph belongs to is determined by the r-cuts 
which may be made on it, as follows: 

C\. No r-cut may be made on the diagram, i.e. it is r-particle irreducible; 

C2: At least one internal r-cut may be made on the diagram and no mixed or final 
r-cut may be made. 

C 3 : At least one initial r-cut is possible on the diagram, but all other types of cut 
are not possible; 

C4: At least one mixed r-cut is possible, but a final r-cut is impossible; 
C5: At least one final r-cut may be made. 

This classification is motivated by where the "latest" r-cut may be made on a particular 
graph. If it is an internal cut then the graph must be placed in C2, but if it is an external 
cut then the graph is placed in one of C3-C5 depending on what sort of cut that "latest" 
r-cut may be. 
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Because each graph contributing to A belongs to one and only one of these classes it 
follows that A may be expressed as the sum of the five expressions found by summing the 
individual classes separately. 

The sum of class C\ is clearly the connected s-channel r-particle irreducible m — > n 
amplitude, A^)_ m . Each of the classes C 2 to C 5 may be summed using the following result, 
which is known as the last internal cut lemma (LICL): 

Lemma 1 (Last Internal Cut) Any (r — l)-particle irreducible diagram which admits an 
internal r-cut has a unique internal r-cut which is nearest to the final state. 

For a proof of this result see paper I. Explicit results for the sum of classes C3-C5 are also 
given there. Note also the following results, more fully explained in paper I: 

1. The Last Internal Cut Lemma (LICL) may not be applied directly when constructing 
the sum of classes C3-C5; 

2. In order to correctly sum classes C 3 , C 4 and C 5 one not only needs to restrict the 
s-channel cut-structure of the amplitudes in the sums of these classes, but also to 
place constraints upon the cut structure of these amplitudes in channels other than 
the s-channel. (This is a point whose full implications were apparently not realized by 
Taylor himself.) However, many of these additional restrictions prove to be trivially 
satisfied by the amplitudes under consideration, due to the stringent nature of the 
s-channel irreducibilities enforced by the LICL. 

B. The double-counting problem and its solution 

As was also explained in paper I, the adoption of Taylor's procedure for summing the 
classes C3 to C5 leads to the overcounting of certain diagrams, as follows. 

In the case where n < r there are complications with the use of the last internal cut 
lemma in classes C 4 and C 5 . If m < r similar difficulties arise in classes C 3 and C 4 . Taylor's 
solution to these difficulties is to split any class C in which such a problem occurs into 
sub-classes S. These sub-classes are defined in such a way that the last internal cut lemma 
may be applied to the diagrams within them without complications, and consequently their 
sum may be constructed correctly. But, Taylor attempts to then construct the sum of C by 
summing over all possible sub-classes S. This is correct, provided that the sub-classes are 
defined so as to be completely disjoint. However, with Taylor's definition of sub-classes of 
C the sub-classes are not disjoint. Certain diagrams belong to more than one sub-class and 
so are double-counted. In paper I the details of this problem were discussed and specific 
examples given. It was then shown how to avoid the double-counting problem by using the 
following procedure to correctly construct the sum of any class C in which double-counting 
is a possibility. 

The procedure involves considering each sub-class S of C in turn. The sub-classes are 
added to a "running tally" of C, one by one. However, a particular sub-class S may only 
be added to the running tally after all of the diagrams within S which have already been 
included in the "running tally" have been eliminated from S. 

We express this mathematically by saying that if C is the union of n sub-classes, i.e. 
C = U'1 =1 Si, then the correct sum of C is: 
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n n n 

J2 <sins , i + E SinSjHSk-.... (io) 

i=1 *, J = 1 i,j,k = l 

i < j i < j < k 

The difficult part of the application of this formula is the discovery of which diagrams 
are in the intersections Si fl Sj, Si R Sj fl Sk, etc.. Because of the way the sub-classes are 
defined any diagram which is in a set Si fl Sj must admit two possible "unique" latest r-cuts, 
each cutting a different set of external lines. Consequently, the diagrammatic content of the 
set Si fl Sj may be found by examining the full s-channel cut-structure of the amplitude that 
is the sum of S% and so determining which portions of that amplitude have a cut-structure 
which means that they also belong to sub-class Sj. In covariant perturbation theory, or, 
for that matter, in any time-dependent perturbation theory, the result of this examination 
will depend, not only on the s-channel cut-structure of the sub-amplitudes making up the 
amplitude representing Si, but also on their cut-structure in other channels. Therefore, 
once the double-counting is eliminated by doing the subtraction recommended in Eq. flTU| ) 
certain constraints are effectively placed on the cut-structure of these sub-amplitudes in 
channels other than the s-channel. Thus, in spite of the fact that we consider only the 
s-channel cut-structure of the full amplitude, the elimination of double- counting effectively 
places constraints on the cut-structure of the sub-amplitudes in channels other than the 
s-channel. 

This modification to the Taylor method provides a way of deriving covariant four- 
dimensional scattering equations which are free of double-counting for any system in which 
it is possible to obtain a graphical perturbation expansion for the relevant amplitudes. The 
only caveat is that the particles involved must be fully dressed before the Taylor method is 
applied. In Appendix |A] we explain how this is accomplished for a system of nucleons and 
pions. Also, in Appendix [B]the Taylor method is used to derive equations for the amplitudes 
which are input to the NN — -kNN equations: the tcN amplitude and the wNN vertex. In 
the next section we demonstrate how the Taylor method as outlined here would be applied 
to the NN t-matrix, Tjvjv- 



IV. THE N — N AMPLITUDE, PART 1 

Now consider the N — N amplitude, Taw- In this work we are interested in particle- 
particle scattering, and so we may consider the particles in the initial and final states to 
both be nucleons, not anti-nucleons. Consequently, the baryon number in any intermediate 
state must be two, and so the amplitude T/vv can only have an s-channel one-particle 
irreducible part T^ N . Furthermore, since, by assumption, both nucleons are fully-dressed, 
this amplitude must be connected. Taylor's method may therefore be applied to this T^ N 
in order to derive an integral equation for NN scattering. 

So consider any diagram contributing to T^ N . Any cuts which could be made to place 
this diagram in C3-C5 would merely expose dressing contributions on the external nucleon 
legs. Since all particles are fully dressed such dressing contributions must be zero, and hence 
classes C3-C5 are empty. Meanwhile, class C\ sums to the connected two-particle irreducible 
NN amplitude T$ N , and class C2 yields, using the LICL, the product T^ N d\d2T^ N where d\ 
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and g?2 are the fully-dressed free propagators for nucleons 1 and 2. Hence, Taylor's method 
gives: 

T (l) _ T (2) T (2) , , T (l) (u s 

Eq. (|Tl|) is, of course, a Bethe-Salpeter equation for Tjp N , with T^ N the kernel of the 
equation, containing the diagrams to be iterated in order to produce the full amplitude. 

(2) 

However, above the pion-production threshold Tjyjy contains part of the effects of inelasticity. 
(The fully-dressed nucleon propagators also contain some inelasticity effects.) Consequently, 
if the NN sector is to be modelled correctly it is necessary that we examine the NNir cut 

(2) (2) 

structure of T^jy. However, if this three-particle cut structure of T N ^ is to be examined via 
the Taylor method we have m = n = 2 and r = 3, so the conditions for double- counting in 
classes C 3 , C 4 and C 5 are satisfied. Therefore, we must proceed with caution. 

(2) 

Applying Taylor's method directly to T N h means that we sum each class individually, 
as follows: 

C\. The sum of class C\ is clearly the connected 3PI 2^2 amplitude, T$ N . 
C^- Using the last-internal cut lemma we may show that the sum of class Ci is: 

F^dM^F^, (12) 

where F^) is the connected r-particle irreducible NNtt — > NN amplitude, and 
F^ is the connected r-particle irreducible NN — > NNtt amplitude. Now, C 2 
may not contain any diagrams admitting a three-cut which involves a final-state 
line. It might be thought that in order to enforce this condition in Eq. (|T2]) 
certain restrictions would need to be placed on F^ 3 '. However, when these extra 
conditions are examined closely it is found that they are all satisfied trivially, 
due to the s-channel three-particle irreducibility of F^ 3 \ (This is, in fact, a 
special case of a more general result, that in C 2 such extra conditions are always 
unnecessary, because the cuts which may intersect a final-state line and so place 
diagrams in C4 or C5 are precluded by the s-channel cut-structure imposed on 
the amplitudes involved.) 

C3: As explained in paper I, the sum of class C3 must be constructed with care. This 
is because the LICL applies only to internal cuts, and so cannot be used directly 
on diagrams in C 3 , since, by definition, diagrams in C3 may have cuts which 
intersect both internal and external lines. 

In order to overcome this problem Taylor constructs two sub-classes of C3, one 
of which contains those diagrams in which the "last" cut cuts the line Nl, and 
the other of which contains those diagrams in which the "last" cut cuts the line 
N2. (Here the sense in which the cut we are talking about is the "last" cut was 
defined in paper I.) Taylor calls the first sub-class C^ N1 ^ and the second C^ N2 \ 
Using the LICL he constructs their sums: 

d m} = F®d 2 d n fU\l), (13) 
Cf 2} = F^d 1 dJ^(2), (14) 
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where (i) is the 1PI nNN vertex corresponding to the emission of a pion by 
nucleon i while nucleon i is a spectator. I.e., i is defined by: 




(15) 



(Note that once again it might be expected that conditions must be placed on 
_p(3) m orc l er to ensure that no diagram which actually belongs in C4 or C5 is 
included in C3. And, indeed, as we shall see below, in some cases the amplitudes 
used in C3 do have to be restricted. But, in the case under discussion here 
the s-channel 3PI of guarantees that only diagrams which belong in C3 are 
included in the sums fll3l) and fll4]).) 

Taylor then claims that the sum of C3 is merely the sum of the two sub-class 
sums (pT3|) and (0). But this is false, since certain diagrams will have two distinct 
possible "latest" cuts, one cutting Nl and one cutting N2. (See, e.g. Figure 0.) 
These diagrams will belong to both sub-classes, and so will be double-counted 
in such a summation. In order to correct this double-counting we apply Eq. (|IT|) 
so as to obtain the corrected sum of class C 3 : 

C 3 = j2F^d;d n fV\?)-D 3 , (16) 

i=l 

where: 

D 3 = cf 1} ncf 2} , (17) 



has been subtracted in order to remove the double-counting. 

C4: Adopting a similar approach to that used to sum C3, we follow Taylor, and divide 
C4 into two sub-classes: C\ N2 ^ N1 \ which contains all those diagrams in which 
the "latest" cut cuts lines N2' and Nl, and C\ Nl ^ iV2 ^ j which contains all those 
diagrams in which the "latest" cut cuts lines Nl' and N2. The last internal cut 
lemma may then be employed in order to write: 

Cr }{m} = fV(2)dJ^(l), (18) 
Cl N1 ' }{m} = fW(l)dJU\2). (19) 



Note that in order to stop diagrams which should be in C5 straying into C4 we 
have had to restrict the amplitude f^ 1 ^ used in ( |T8"D and ( |T9"D to only contain 
diagrams which are 2PI in the N' <— N + 7r'-channel. We have denoted the 
resulting amplitude by f^. If the restriction on is not enforced then 
the three-cuts depicted in Figure ^| may be made on certain diagrams in C4, 
indicating that these diagrams actually belong in C5, not C4 

If we now apply Eq. (0) we find that the correct sum of C4 is: 
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c 4 = E /%x/ (1)t (*)4- - A, 



(20) 



with: 

d 4 = cl N2 ' }{m} n cl m ' }{N2} . (21) 

C5: Once more this class is split into sub-classes. This time the sub-classes are C^ N1 ^ 
and C5 , which contain, respectively, all diagrams in which the "latest" cut 
cuts lines NV and NT!. The LICL gives: 

C { 5 Nj ' } = f^\j)d 3 d^\ (22) 
for j = 1,2. Consequently we have: 

^5 = E/ (2) (j)^ (2)t -^5, (23) 

3=1 

with: 

D 5 = cf 1,} ncf 2 ' } . (24) 

If we combine the above results for the sums of the individual classes we find: 
T£ ] N = T®> + F^d 1 d 2 d 7r (F^ + E f {1) \W) ~ Ds 

i 

+ E/ (2) 0')^^(E^/ (1)t (^)< rl + ^ (2)t ) " Ai - D 5 (25) 

j i 

where here, and throughout the rest of the paper, the sums over i and j are understood to 
run over i,j = 1,2. This equation, but without the restriction on in class C4 and the 
subtractions for double-counting, was also derived by Avishai and Mizutani [ffO] and Afnan 



and Blankleider ||1T]J17| . Below we shall calculate the precise values of the subtraction terms 



(2) 

D3, D4 and D5. But, in order to achieve this and so derive an equation for T^y we will 
need to know the structure of the amplitudes in Eq. (|25|). 

Firstly, consider the amplitudes T^ N and . Explicit examination of the structure 
of these amplitudes would involve exposing NN7m intermediate states. Since our primary 
concern here is the derivation of equations which treat the NNw intermediate states of the 
theory correctly, we argue that these amplitudes may be safely neglected or approximated. 
The equations resulting from this approximation will then satisfy iViV and NNir unitarity, 
but not NNnir unitarity. As a first approximation we set F( 3 ) and T^ N to zero, in order to 
make our equations as simple as possible. This approximation is, in fact, exact in the case 
of an Nn interaction Lagrangian containing only a nNN vertex. 

The use of this approximation means that the sum of both sub-classes of C3 is zero, and 
so D 3 = 0. Therefore, given this assumption, T^ N obeys: 
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= Y,f (2) U)d-AF^ + J25 ij f®(j)dJM(i) - D 4 - D 5 , (26) 

j ij 

with D4 and D5 given by equations (|21| ) and fl2~4|) respectively. The amplitude J*- 2 ** which 
appears in (^) can be extracted from the model Lagrangian, as demonstrated in Appendix 
[B|. As far as the amplitude /Wt is concerned, for the moment we merely observe that it 
must contain some subset of the diagrams which are summed by solving Eq. ( |B2"D for 
Therefore the only amplitude whose structure remains to be investigated is . 



V. THE TWO-PARTICLE IRREDUCIBLE, NN -> NNtt AMPLITUDE, 

(2) 

Our examination of the two-particle irreducible NN amplitude T^jy- revealed that the 
equation it obeyed involved the connected 2PI iViV — ► NNtt amplitude F^ . In this section 
we examine this amplitude using the original Taylor method, and eliminate the double- 
counting inherent in that method, using the modification to the Taylor method described in 



Section 1IIB above. 



In this case the sum of the Taylor classes and sub-classes is as follows: 





= F (3)t 


(27) 


c 2 


= M^dxd^F^ 


(28) 


c {Nl} 


= M®d 2 djW\l) 


(29) 


C {N2} 


= M 2 (3) d 1 rf 7r /« t (2) 


(30) 


( j{Nl'}{N2} 


= fi(l)^/ (1)t (2) 


(31) 


^{JV2'}{JV1} 


= f&(2)d ir /« t (l) 


(32) 


C W}{N1} 


= T$Uf {1) \l) 


(33) 


C W}{N2} 


= r£W (1)t (2) 


(34) 


C {Nl'} 


= t^ ) N (l)d 2 d 7T F^ 


(35) 


C {N2'} 


= tW(2)did«F^ 


(36) 




= T$Ud 2 FW\ 


(37) 



where we have labeled the incoming nucleons Nl and N2, the outgoing nucleons NV and 
N2', the outgoing pion tt', and t^(i) is the 2PI tt — N t-matrix with nucleon i as a spectator. 

Here, and M2 are all connected 3PI NNtt NNir amplitudes. However, 

(3) 

constraints beyond simple s-channel 3PI have been placed on the amplitudes M{ and 

(3) (3) (3) 

M2 . These constraints mean that M[ ' and M£' are different and neither is equal to the 
sum of all s-channel 3PI NNtt — > NNtt diagrams, M^ 3 \ It is found that in order to stop 
diagrams actually belonging to C 4 being mistakenly included in C 3 we must define M- 3 ^ to 
be 1PI in both the channels: 



Ni' + Ni «- Ni' + Ni + tt' + tt 
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and: 

7r' + Ni <- Ni + Ni + Ni' + 7r. 

It might appear that conditions beyond s-channel 3PI would also have to be placed on the 
appearing in C 2} in order to prevent diagrams from C4 and C5 being included there too. 
However, in the case under consideration in this section the possibility of any cut placing 
a diagram which is summed in C 2 above in C4 or C5 is precluded by the s-channel 3PI of 
Therefore, as was the case with the amplitude F^ 3 ^ used in C 2 and C 3 in the previous 
section, no constraints beyond s-channel 3PI are necessary. 

Note that once these restrictions are imposed on the problem discussed by 

Kvinikhidze and Blankleider [EH], where diagrams hidden in M®did n f m \i) actually rep- 



resent processes included in C4, does not occur, since the diagrams in C4 which give rise to 
these difficulties are specifically excluded from C^ N1 ^ and c\ N2 ^ by the above irreducibil- 
ity conditions. Therefore we may safely conclude that all diagrams in C 2 and C3 will only 
represent processes involving some piece of the connected s-channel 3PI NNtc — > NNn am- 
plitude. In this paper any such process is said to involve a "pure three-body force" . In the 
calculation here we assume that to a first approximation these pure three-body forces are 
negligible, by an appeal to the same argument which resulted in us neglecting the amplitudes 
and above. Consequently, we assume that the sum of classes C 2 and C3 is zero. 

The inclusion of a three-body force, in our calculation, would result in considerable 
changes to the ensuing argument, since at numerous stages during our discussion below we 
have appealed to the absence of a three-body force as a reason for eliminating certain double- 
counting subtraction terms. Nevertheless, integral equations for the amplitudes involved may 
still be derived, even if ^ 0. In that case, however, it is considerably harder to see how 
such equations may be manipulated into a form for practical calculation. 

Once these observations have been made it is a simple matter to sum the above expres- 
sions for the sub-classes in order to obtain the following equation for : 

F^ = ^t {2 \a)d- 1 d l d 2 d n (J2Siaf {1) \^ + F^) - D A (38) 

with the sum over a running over a = 1,2,3 here and throughout the rest of this work; 
(a) being defined by: 

t«(a) = /*^ if« = i = l,2 (39) 
[ Tfrh if a = 3. 

and: 

D A = c 4 {7r ' }{Afl} n d n ' }{N2} + d N2 ' }{m} n cl n ' }{N2} + cl m ' }{m} n cj 7r ' }{iV1} 

+ C {N2'}{N1} n C W}{N1} + C {N1'}{N2} R c {n'}{N2} + C {N1'}{N2} R C {N2'}{N1} ^ 

included in order to eliminate the double-counting in C 4 , as is prescribed in Eq. ([101) . Note 
three things about this result: 

1. The D4 here is obviously different to the D 4 defined by Eq. fl2"Ip in the previous section, 
just as the sums of the Taylor classes in this section are different since we are dealing 
with a different amplitude here. 
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2. The maximum number of final-state lines cut by any three-cut is Sf = 1. Hence the 
condition n > 2s t holds and this condition is sufficient to guarantee the absence of 
double-counting in C 5 . (See paper I for details.) 

3. All three-sub-class intersections are zero in this particular case, thus is given by 
the sum of the two-sub-class intersections listed here. 

We now calculate the six two-set intersections involved in D^. In order to do this we note 
that any diagram which is in both C\ and C\ must admit two "latest" three-cuts, 
the first of which cuts the three lines A', B and one other, and the second of which cuts the 
lines E , F and one other. Consequently, the intersection of C\ and C\ may be 
found by investigating which diagrams contributing to the sum C\ admit a three-cut 

involving E', F and one other particle. 

Looking at the sum of ct' }{N1} we see that no three-cut involving it', N2 and one other 
particle is possible. Therefore we conclude: 

cf }{m} n cf }{N2} = o. (4i) 

Similarly: 

cf 2 ' }{m} n cP {N1} = 0, (42) 

and: 

c {m>} { N2} n c w H N2 } = Q (43) 

Now consider c\ N2 ^ N1 \ Upon examining the diagrammatic sum of C\ N2 ^ N1 ^ we ob- 
serve that there are a number of diagrams contributing to the sum in which a three-cut 
involving n' and N2 may be made (see Fig. ^). Such a cut may be made in the portion of 
C\ N2 ^ N1 ^ which sums to: 

t&^K/WV), (44) 

where v* N is the u-channel 1PR part of t^. In Appendix it is shown that the part of 
t£jv which is one-particle reducible in the s-channel is the pole diagram, d^f^- By 
changing the argument of Appendix ^ to apply to the it-channel instead of the s-channel it 
may be shown that v* N is merely the crossed term: 

«& = / (1) ^/ (1)t - (45) 

Therefore, 

C {N2> }{ m } n Q K }{J V2} = v X (2)rf7r/ (l)t (1) (46) 

= ^OP^i/ (1)t (2), (47) 
Note that here fops is the full one-pion exchange potential: 
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v OP E = f m (2)d 7r fM\l). (48) 

This argument makes it explicit how the term (H3) comes to be in two sub-classes. It is this 
inclusion of (|46|) in two different sub-classes which leads to the double-counting pointed out 
by Kowalski et al. ^ . 

Naturally exactly the same argument, but with the roles of iVl and N2 reversed, may 
be used to show that: 

C {NV }{ N2 } n C WHN1 } = v X {l)dwf (^ {2) (4Q) 

= v PEd 2 f {l) \l). (50) 

The identification in the last line may be made because the one-pion exchange potential 
may be written as: 

v PE = f (1) (l)d^ 1) \2), (51) 
since we are working in a time-dependent perturbation theory, and so: 

/Wt = / (i). (52) 



Note that in order for Eq. (p2|) to hold all particles involved must be fully dressed, so that 
is one-particle irreducible in all channels. 

Finally we attempt to construct the intersection of C\ Nl ^ N2 ^ and C\ N2 ^ N1 ^ . Any 
diagram which is in C\ N1 ^ N2 ^ and admits an s-channel three-cut involving N2', Nl and 
one other particle will also be in C\ m ^ N1 \ As is shown in Fig. || such a cut is possible if 
the amplitude t^' N used in C\ N1 is one-particle reducible in the i-channel. Therefore in 
models including diagrams such as Fig. || the intersection will be given by: 

cr }{N2} n cf 2 ' }{m} = <^(i)4/ (1)t (2), (53) 

where v% N is the t-channel one-particle reducible part of the s-channel 2PI 7r-N amplitude. 
However, in our model pions do not interact with any particle other than nucleons, and so 
the n — N t-matrix is automatically one-particle irreducible in the t-channel. Consequently, 
the cut shown in Fig. || is simply not possible, leading to: 

c{ m> H N2 } R C {m> }{ m } = Q (54) 

Therefore we have: 

Da = E / (1) (0d// (1)t (04r/ (1)t 0')^ (55) 

ij 

= J2v x (i)d„fM\j)5 ij (56) 

ij 

= Y. v oPEd 3 f 1] \i)6 l3 (57) 



'j 



As these rewritings show, when this result for D 4 is substituted into the above integral 
equation for the subtracted diagrams may either be used to modify the n — N t-matrix 
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or the N — N t-matrix. As will be seen below, which t-matrix we choose to modify here will 
affect our later results for other amplitudes. In this work we choose the former approach, 
however we stress that modifying the N — N t-matrix is equally legitimate. When these 
subtractions are made we find that: 

F (2)t = j2 v R (a)d- l 5 m d^f {1) \i) + t(2) {a)d- 1 d 1 d 2 d„F^ (58) 

ai ol 

with: 



v R (a) = \ t «$>- V ® *° = i =1 ' 2 (59) 



, T^n if a = 3. 

Note that: 

1. If we choose to modify the NN t-matrix Eq. (|58"D still holds, but with v R (a) given by: 

v R (a) = l t( f 2 ^ if « = z = 1,2 (6Q) 

[ T^jv - vqpe if ol = 3. 



(2) (2) 

2. Provided we use amplitudes T^n and which themselves contain no double-counting 
the integral equation for , Eq. (|58|), is completely correct and, in turn, contains 
no double-counting. 

(2) 

3. There is a bootstrap problem here. We set out to determine T^jv and found it depended 
on F^\ Now we find that F^ ] depends on T$ N \ 

However, leaving this bootstrap problem aside for the moment, Eq. ( j58|) may be iterated 
in order to sum the multiple-scattering series for F^ . If this is done and the definitions 
fl39f)-(l65D used, we obtain: 

F (2)t = Y,f 1 \(3)dp 1 d l d 2 dJp l fV\ i )d; 1 

p% 

+ Y,t [1 \a)d- 1 d 1 d 2 dM$d 1 d 2 dJ {1 \p)dp 1 d l d 2 dJp^ (61) 

a.p% 



where the U^Js obey the four-dimensional covariant AGS equations plj : 

U«p = Sapd^d^d- 1 + VW^di^rf^?, (62) 

7 

= Sapd^d^d- 1 + J2 U^d^d^i^d- 1 ^, (63) 

7 

and: 

€n = (1 + #i«>£v = v R N + t^ N d N djH (64) 
t^NN — TpfN, (65) 
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where T^ N is, in principle, the full NN t-matrix, which is related to T^ N by Eq. (ITU) . 

Note that all two-body amplitudes in this series are the full s-channel 1PI two-body 
amplitudes except that the first tc — N scattering after the pion emission involves a 
t-matrix, t^y, which has had the crossed term removed from the potential. This removal of 
the crossed term is the only way in which this result differs from that obtained by Avishai 
and Mizutani and Afnan and Blankleider. It is necessary because, as was first pointed out 



by Kowalski et al. P0 |, the inclusion of the crossed diagram in the first n — N amplitude after 
pion emission leads to double-counting of terms such as Eq. (f46[) and (|49|) if the calculation 
is done in a time-dependent perturbation theory. Because we have used the modified Taylor 
method, which is specifically designed not to lead to double-counting, the above equation 
does not contain such double-counting, even though it was derived in a time-dependent 
perturbation theory. 



VI. THE N-N AMPLITUDE, PART 2 

Having obtained the correct expression for F^ 2 ^ we may now return to our calculation 
of Ttflf. Recall that above we derived Eq. ( PB[ ) for T^ N , with the first term generated by C 5 
and the second by C 4 . However, recall also that both terms were expected to contain double- 
counting. Having derived an expression for F^ which is itself free of double-counting we 
may now determine the factors _D 4 and D 5 , which were introduced into Eq. (|26|) in order to 
compensate for this double-counting. 



A. Calculating 

Firstly, note that: 

d, = c\ NV]{m] n cl N2 ' }{N1} , (66) 

with C\ Nl ^ N2 ^ and C\ N2 given by the expressions (0) and (|19f) . 

Now observe that since all particles are fully dressed all vertices are at least 1PI in all 
channels. Therefore in a time-dependent perturbation theory: 

/(i)t = / ( 2 ). (67) 

Consequently: 

f^(2)dj^(l) = f^{l)dj^{2), (68) 

i.e. one sub-class merely reproduces the sum of the other. Topologically this happens because 
the two cuts N1'ttjN2 and N2 f 7TjNl (here / stands for intermediate line) may be made on 
any diagram in C4 and so all diagrams in C4 belong to both sub-classes. Therefore: 

D A = cl NV}{N2} = /< a >(2K/M(l) = C[ N2}{N1 ' } = f^(l)dJW(2). (69) 
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B. Calculating D§ 
In Section [TV] we denned D 5 to be: 

d 5 = c { 5 N1 ' } n d m ' } . (70) 



In order to discover which diagrams are actually in D 5 we examine the sum of a 



{m>}. 

5 



Ci m ' } = f^\l)d 2 d^\ (71) 

and note which diagrams in c\ Nl ^ admit a three-cut involving N2' and so should be included 
in D$. Naturally, provided that we use an expression for which is itself free of double- 
counting, such as Eq. fl6"T|) derived above, the sum of C^ 1 * will itself be free of double- 
counting. Consequently, such a procedure will correctly identify all of the diagrams which 
must be included for subtraction in D 5 . 

Now substituting in ( |7lD for F^ from Eq. (|6l[) and iterating using the AGS equations 
gives: 

Cf V} = E^^lWtf'WCW-/ 1111 ^) + f {2) (l)d 2 dJ^ N (l)d n f^\2) 

ai 

+ ^/^(l^dA^cO^i^ (72) 
We write: 

i{Nl>] 



CI = ci + c 2 + c 2 + c 3 + c 4 + c 5 + c 6 , (73) 



where: 



c 1 = /«(l)djS(2)d 2 ^/( 1 ) t (l), (74) 

c 2 = / (1) (l)^W (1)t (l)> (75) 

c 3 = / (1) (l)rf2^ri 1 irf 1 /« t (2), (76) 

c 4 = / (2) (l)^S(l)^/ (1)t (2), (77) 

c 5 = E / (1 Hl)^^tS(2)rfi^^g ) rf 1 rf 2 4t (1 H/5)^ 1 ^^/ (1)t (^), (78) 

pi 

c 6 = E / (1 Hl)rf2^4 1 ^i4f/g ) rfirf 2 ^t (1) (/5)^ 1 ^^/ (1) ^). (79) 
It is clear that: 

D 5 = £c n C 5 {iV2 ' } . (80) 

a=l 

Therefore in order to determine the value of D 5 it is necessary to examine the six terms 
Ci-Cg individually, probing carefully in order to determine which diagrams from each term 
admit a three-cut involving N2' . 
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1. Cl 

We wish to calculate: 

Cl ncf 2 '>. (81) 

In order to do this we observe that if the diagram representing c\ is to admit an s-channel 
three-cut involving N2' and two internal lines then the amplitude t^h will have to be two- 
particle reducible in the N' <— tt' + tt + N channel (see Fig. |7|). Therefore we investigate t^ N 
in order to discover what portion of it is two-particle reducible in this channel. We begin by 
using Eq. ( |S4D to write: 

t£& = v R + v x d N dJ^ N + v R d N dJ^ ] N . (82) 

The second term here is 2PR in the channel N' <— tt' + tt + N. Therefore, substituting this 
term into Eq. ([T4]) for c\ it is found that the diagram: 

/«(l)^/i 1 )(2)rf 1 /« t (2)d 1 d4 1 i(2)d 2 d 7r /« t (l) (83) 

is in ci fl c\ N2 \ 

Examining (j32"D it is now seen that the only other way in which tfy could be 2PR in the 
N' <— tt' + 7r + A^-channel is if v R is 2PR in the same channel. Therefore we now examine 
the effect of placing the N' <— n' + tt + A^-channel 2PR part of v R in the expression for c±. 

Claim 1 The portions of the expression: 

f 1 \l)d n v^(2)(l + d 1 dJ^ N (2))d 2 d n f^\l), 

where is the u-channel 1PI, s-channel 2PI and N' <— tt' + tt + N -channel 2PR tt — N 
interaction, which are also in Cl , given the assumptions of this calculation, are: 

/SHl)^J^H2)rfi^(2)^ 1 rfi/« t (2)^ 2 4/«V). (84) 

and: 

+/gHl)^*/SH2)^di*£U^/S )t (2)d«/ 1 g )t (l) (85) 



For a proof of this result see Appendix |C[ 

We never intended to include the possibility of tt— tt interaction in our equations, since we 
think it is too small to make a significant difference to the calculation. This is a reasonable 
approximation, since in a model with only pions and nucleons two pions may only interact 
via the exchange of nucleon-anti-nucleon pairs. Consequently the lowest-order diagrams for 
tt — tt scattering in a model with a ttNN vertex and a tt — N contact term are the two 
diagrams shown in Figure When these diagrams are used in (|8~5|) they produce a NNNN 
intermediate state. Considering that we are already neglecting explicit ttttNN states in 
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the calculation, and that we are also neglecting processes involving mesons other than the 
pion, e.g. the rho meson, the error made in ignoring diagrams such as those represented by 
Eq. ( pop is insignificant, and consequently from now on we neglect the expression (R3J) as a 
possible source of double-counting. 

Therefore the sum of all diagrams in c\ fl C\ is: 

/i! ) (l)^/( 1 )(2)d 1 tf )t (2)rf 1 ^ 1 i(2)rf 2 4/i 1)t (l) 
+/SHl)^jSH2)rfitS,v(2)rf 1 <^) t (2)^ 2 rf 2 /£) t (l). (86) 



2. c 2 

The term c 2 has one possible cut involving N2' and two internal lines, as shown in 
Figure |9|. This diagram shows that any portion of T^jy which is two-particle reducible in the 
NV <— iVl + N2 + iV2'-channel will, when placed in the expression for c 2 , lead to diagrams 
which are in c 2 fl C^ N2 \ 

Consequently, we need to identify all NV <— iVl + N2 + iV2'-channel 2PR portions of 



Txlf. In order to do this write: 

(1) _ ™(2) . ^(2) , , ^(1) 



Tnn — ^nn + T NN d\d2T NN . (87) 

T^jy is then split into two pieces: the t-channel 1PI part, which is denoted by T^ N , and the 
t-channel one-particle reducible part, which is clearly just the one-pion exchange potential. 
Therefore we obtain: 

Tnn — v ope + T^n, (88) 
(note that this is exactly the alternative decomposition given above in Eq. (j60|) ) leading to: 

Tip N = /«(2)d 7r /« t (l) + Tjfl + f {1 \2)d^f^\l)d 1 d 2 T^ ) N + f^d.d^l (89) 

Examination of the first term on the right-hand side of this equation shows that if f^(2) 
is two-particle reducible in the N' <— iW channel then this part of T$ N will be 2PR in the 
NV <- Nl + N2 + 7V2'-channel. Use of the LICL then shows that: 

/^(l)d 2 ^/i 2 )(2)d 1 4^(2)4/« t (l)d 2 /« t (l), (90) 

is in c 2 PI C^ N2 ^ and hence in D 5 . 

Furthermore, the third term on the right-hand side of (|89"D is completely 2PR in the 
NV <— Nl + N2 + iV2'-channel and thus the diagram: 

/W(l)^/i 1) (2)^/i 1) V)^^Wi 1)t (l). (9 1 ) 
is in the set c 2 fl C\ N2 ^ . 
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Now consider the second and fourth terms of the right-hand side of Eq. fl39|). The only 
way in which these terms lead to portions of T^ N which are 2PR in the Nl' <— N1+N2+N2'- 
channel is if T^ N is itself 2PR in the same channel. Arguments given in Appendix |B| show 
that the only contributions of such diagrams to c 2 fl C5 are: 

/i^2)4j£)(l)d 2 ^(l)d 2 4 2 /S )t (l)rfi^jiI )t (2), (92) 

and: 

/(;)(2)d 1 </«(l)^ 2 d 2 T^d 2 ^) t (l)d 1 <^) t (2), (93) 

where T^lf is not only one-particle irreducible in the s-channel, but also 1PI in the t-channel. 

Therefore, the parts of c 2 which are also in C§ m ^ and so are included in .D5 are the 
diagrams corresponding to Eqs. Q9"0|), (|9~T|), ( |9"2"| ) and (|9~3"D. 



3. c 3 

Examination of the diagram representing this term shows that two three-cuts involving 
N2' are possible, as shown in Figure |H]. The possibility of these three-cuts means that any 
diagram contributing to C3 which contains a part of the NN t-matrix T^ N which is: 

1. 1PR in the ^-channel; 

2. 1PI in the t-channel but 2PR in the Nl' <- Nl + N2 + iV2'-channel, 
is in C3 fl C^ N2 \ 

As far as possibility 1 goes, substituting the 1PR piece of T^ N , vope-, into the expression 
([TB|) for C3 indicates that the term: 

/SH^^/ff^^/ffVjdi/ff^ 2 ). ( 94 ) 

is in c 3 nCf' } C D 5 . 

We also showed above and in Appendix [B] that the Nl' <— Nl + N2 + iV2'-channel 2PR 
parts of the t-channel 1PI NN t-matrix are: 

1. fM(2)d n fW\l)d 1 d a T$} f] 

2. f^{2)d 1 dj^ ] 

3. f^(2)d 1 d n F^d 1 d 2 T^ N ; 

where the amplitude is defined in Appendix [D|. However, when these three diagrams 
are substituted into the expression for C3 only diagram 1 and the portion: 

f^(2)d 1 d^ N d 2 f^\l) (95) 
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of diagram 2 produce diagrams which are in c 3 fl CI . This is because all the other 
diagrams produced by the substitution of NV *- Nl + N2 + iV2'-channel 2PR parts of f^ N 
into C3 will only be contained in C5 if a three-body force is included in the calculation. 
Therefore the contents of C3 fl C^ N2 ^ are the diagram (0), and the diagrams 

/«(l)^ 2 /^(2)^tf )t (l)d 1 d 2 Ti 1 i ( i 1 /( 1 ) t (2) (96) 

and: 

/£Hl)*/SH2)^TjWd ra cfc/W t (l)'ii/S )t ( 2 )- ( 97 ) 



4. c 4 



If the three-cut drawn in Figure [TT| may be made then the diagram representing C4 will 



{N2 1 } 



also be in C L 

Therefore diagrams which sum to: 



f^{l)d 2 d^{l)d^ R {2) 



(i) t. 



(98) 



where f^pR (2) is the N' <— iV-|-7r-channel 2PR part of the IPI tcNN vertex, are in C4<~)C 5 
and so are in D 5 . 



{N2'} 



5. c 5 

Examination of the diagram representing the sum of C5 shows that an s-channel three-cut 



involving N2' cannot be produced by cutting lines beyond the first interaction (see Fig. |12| ) 

(2) (2) 
Therefore, the structure of and the amplitudes appearing after are irrelevant to 



the calculation of C5 fl C^ m \ 

Now an s-channel three-cut involving N2' will be possible on this diagram if: 

1. The n — N t-matrix used immediately before the pion absorption is u-channel 
IPR; 

2. The w-channel IPI part of that tfy is N' <- N + vr + yr'-channel 2PR. 

Using the above results for the -u-channel IPR part of and the N' <— N + it + n'- 
channel 2PR part of the u-channel IPI t^, but ignoring those diagrams which are only in 
CI if a three-body force is included in the calculation, shows that: 



c 5 n Cf 2 ' } = f$ \i)dA,fU(2)d 1 f£ }t (2)(i + d!^t?i(2)) 

x d 1 d 7r U^d 1 d 2 dJ {1) (f3)d^ 1 5 f3i d- i d n f < - 1) \{}. (99) 
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6. c 6 

As was the case for C5, we observe (see Fig. O) that an s-channel three-cut involving N2' 
cannot be produced by cutting lines beyond the first interaction. Therefore, the structure of 

(2) (2) \N2'\ 

l/y and the amplitudes appearing after l/y are irrelevant to the calculation of c§ fl C\ 
Then an s-channel three-cut involving N2' is possible in Fig. [13] if: 

1. The 1PI N — N t-matrix T^ N used immediately before the pion absorption is t-channel 
1PR. 

2. The t-channel 1PI NN t-matrix, f£l is NV <- Nl + N2 + 7V2'-channel 2PR. 



Using the previous results for these portions of the NN t-matrix and again ignoring those 
:ts of 
we find: 



{N2'\ 

parts of cq which are only also in C5 if a three-body force is included in the calculation 



c 6 n cf 2 ' } = J2 fi 1 Kl)d 2 d^ ) (2)d, / fS ) \l)(l + d 1 d 2 T^ ) N ) 



ft 

x d 1 d 2 U$d 1 d 2 dJ {1 \p)dp% i did 7T f {1) \{). (100) 



7. Overall result for D§ 

Combining the results for the intersections c a fl C5 , a = 1, . . . , 6, as per Eq. flg0"|) 
reveals that: 

+ £/«(l)d 2 4/£ ) (2K^ )t (l)(l + ^^ 

+ f 2 \l)d 2 d^l{l)d^ R f (2) 

+ f$\l)d 2 d 7: ,fP (2)d 1 dJ: ) N (2)d 7T f^\l)d 2 fil ) \l) 
+ B + X + D 1 + D 2 + D n 

+ /( 1 )(l)^ 2 /i! ) (2)^/« t (l)rf lC / 2 T^rf 1 /( 1 ) t (2) 

+ /«(i)M2/5 ) (2)^/^V)^^4 1 id 1 /( i ) t (i) 

+ fi)\l)d n ,fi 1 \2)dj^\2)d 1 d^l(2)d 2 d n f^\l), (101) 
where: 

B = /(: ) (l)rf 2 /W(2) C / 7ri rf 7r2 /« t (l)d 1 /« t (2), (102) 
represents crossed two-pion exchange; X is given by: 

X = /«(l)d 2 /«(2)d 1 ^ i ra^ 2 rf 2 /« t (l)rf 1 /W t (2), (103) 
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the terms Di are: 



A = /«(lK 2 d 2 /«(2)cM^^ (104) 

and is given by: 

Ar = /SHl)^^2^ ) (2)rf 1 ^ 1 T( r 1 irfJ«^2)d 2 /« t (l). (105) 

At this stage we must make a decision. Upon substitution of the expression for D 5 into 
(p6|) many of the terms in D5 can be subtracted from either the N — N t-matrix or the 
n — N t-matrix. In what follows we choose to subtract them from the n — N t-matrix, as 
we did above for the t-matrices in . This has the effect that many of the "potentials" 
generating the ir — N t-matrices, which are, by definition, 2PI in the s-channel, also become 
1PI in the u-channel. In other words, the resultant input ttN t-matrix appearing at some 
points in the equations has the crossed diagram removed from its "potential" . However, we 
stress that it would be equally valid to modify the NN t-matrix. Were this to be done the 
modified NN "potential" would be 2PI in the s-channel and 1PI in the t-channel. That is, 
if we chose to proceed in this way the input NN t-matrix which appears at some points in 
the equations will not contain the one-pion exchange contribution in the driving term of its 
Bethe-Salpeter equation. 

As we shall see further below, the fact that we choose to modify the ti — N t-matrix 

(2) 

means that the ttNN form factor in the one-pion exchange part of T^^ needs to be ad- 
justed. By contrast, if we chose to modify the iViV t-matrix it may be shown that this 
change to the nNN form factor in one-pion exchange is unnecessary. This result is merely 
a specific example of the wider truth that any adjustments which are made in the pion 
production/absorption sector have significant implications for the NN — > NN sector of the 
theory. 



C. Equation for 

Having made this decision we may substitute eqs. (|6lD , (|6^) and ( |101| ) into Eq. (|26"D , and 
simplify using the AGS equations in order to obtain: 

T { ^ N = f 2 \l)dJ^(2) + C d 5 

+ E f {1 \j)S j ad- j dJ- 1 \a)d-%dJ ai f^\z) 

jai 

+ Y,f W U)^d ] dJ^\a)d- 1 d 1 d 2 d^d 1 d 2 dJ^\ 

— Di — D 2 — D n — X — B, (106) 

where: 

4$ = v R + v R d N d^ l \ (107) 

and: 
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£~S = v R + t^d N d w v R , (105 

are both one-particle irreducible in the s and it-channels, but differ in that has the 
crossed term subtracted f 
interaction. By contrast: 

t 1] {t) = i nN = v u + v u d N d n v H + v H d N d 7T t^d N d 7T v H , (109) 



crossed term subtracted from the first interaction, while has it subtracted from the last 

,R , ..Rjj „,R , .JIjj + (i) j. j ..R 



for i — 1,2, has the crossed term subtracted from both the first and the last term. Meanwhile: 

? 1 \a) = Tj}l (no) 

and C| is that portion of the sum of class C5 which contributes to the dressing of one-pion 
exchange and is given by: 

C d 5 =f^(2)d 1 dJ i : ) N (2)d 7T fW\l) 
+f< 2 \l)d 1 d^ ) N (l)dJ^{2) 

-/^(2)d 1 d^W(2)^/^ ) (l)d 2 /« t (l)4^tf )t (l)- (HI) 
Formal manipulations now show that: 

Cl + /( 2 )(1K/Mt(2) = V * pE + DoPEj (112) 

where: 

Dope = fg\l)d n d 2 fg\2)d 1 d^fg\2)d n fff\^ (113) 
^ = / (1) *(i)^/ (1) * t (2); (114) 

with: 

f M* = fW + fWd N d^l (115) 

the vertex which is obtained if the crossed diagram is eliminated (in the manner prescribed 
above) from the % — N t-matrix used to generate the dressed 7rNN vertex. Note that, using 
the definition of t n N, Eq. (|64|), we find: 



^^'-tfWxW (116) 

as depicted in Figure 0. Note that in order to obtain the result (|112|) we have had to use 
the time-dependent perturbation theory identities ( |52"D and (|67|) . 

Resubstituting the result ( |112j ) into Eq. (|106|) yields the final result for t/ 2 -* • 

,(2) 



NN- 

'/;■-; Y orh: -D [ -D 2 -D,-X-D 

l) 



+ E f {1 \j)5 ja d- j djf\a)d- 1 d l dJ ai fW\ 

jai 

+ E f {1 \j)S ja d ] dJ^{a)d- 1 d l d 2 d^d 1 d 2 dJ^ (117) 
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where: 



D^ = D^ + Dope = /i^K^/i^K* d^SU^f (2)d 2 f£\l). (118) 

(2) 

The amplitude Xjyjy is the nucleon-nucleon s-channel two-particle irreducible interaction 
which is to be used in the Bethe-Salpeter equation for NN scattering in this theory, Eq. ( |TTD . 
Hence it plays the role of the nucleon-nucleon "potential" in this work. However, rather 

(2) 

than merely being the sum of a few Feynman diagrams T^jy is the sum of all 2PI Feynman 
diagrams with one explicit pion. Consequently it includes the full explicit coupling to the 
ttNN channel. 

If this result is compared to the equation for T^jy obtained in the derivation of the 
standard NN — ttNN equations [[2}-|T2"|] a number of differences are seen: 

1. The one-pion exchange part of the "potential" has been replaced by the modified 
"potential" : 

V = V* PE -D 1 -D 2 -D n -X-B, (119) 



that actually contains (for reasons to be discussed below) the diagrams D 1 , D 2 , D n , 
X and B which involve multiple pion exchanges. 

In the two terms describing pion rescattering (i.e. the last two terms on the right-hand 
side of Eq. ( |1 1 7| ) ) the tt — N t-matrix has been modified. In the first term the usual 1PI 



~(i) 

n — N t-matrix has been supplanted by the t-matrix t , in which the crossed diagram 
has been removed from both the first and the last pion-nucleon "potential". Likewise, 
in the second term the operators fi 1 ' and govern, respectively, the first/last pion 
rescattering after/before pion emission/absorption. 

All of these changes are the result of our desire to eliminate double-counting from the 
equations. Taking the second set of changes first, in Section |V| it was shown that either the 
ttN or NN amplitude immediately following a pion emission needed to be modified in order 
to avoid the over-counting of diagrams involving the process ([55]). The changes to these 
n — N t-matrices in Eq. (|117|) are made for exactly the same reason. The effect is to place 
constraints on these sub-amplitudes in channels other than the s-channel. As discussed in 
Section (TO], such additional constraints are necessary because in a time- dependent pertur- 
bation theory diagram constraints on the s-channel cut-structure of sub-amplitudes are not 
enough to constrain the overall s-channel cut-structure of the amplitude. 

Turning now to the first set of changes, we divide our discussion into two halves. Firstly, 
let us consider the one-pion exchange piece of V, Vq PE . A very interesting consequence of 
our work here is that the one-pion exchange "potential" Vq PE contains a different vertex 
than that used in the pole part of the it — N t-matrix t^. This is necessary because the 
inclusion of the full ttNN vertex in the one-pion exchange potential leads to double-counting 
with respect to other terms in T^ N . (See Fig. [15] for examples of such diagrams.) These 
other terms are included in the final two terms of Eq. (|117|) . Therefore, they only arise 

(2) 

when the NNtt cut structure of the amplitude Tjyjy is considered. This is a manifestation 
of the fact that the amount of double-counting to be eliminated is entirely dependent upon 
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which unitarity cuts are opened in the analysis. Our above discovery that certain types of 
double-counting do not arise when a three-body force is excluded from the calculation is 
another manifestation of the same effect. Therefore we find that whether or not the vertex 

(2) 

used for the one-pion exchange piece of T^jy needs to be adjusted depends on whether: 
1. the ttNN unitarity cut is opened or not; 



2. the subtraction terms are included in t n N or T NN (see discussion in Section |V1B[). 



Secondly, the other subtractions in V (i.e. D x , D 2 , D n , B and X) are also required 
in order to correct for the overcounting of diagrams in the last two terms of Eq. (|1 1 7|) . 
(Note that once again it is the opening of the nNN cut which has led to double-counting.) 
As is shown in Figs. [16| and |17| Di, D 2 , D^, B and X are all included in this term in 
more than one place. Consequently these five diagrams will be double-counted unless some 
adjustment to the equation for T^ N is made. Because we wish to obtain a set of coupled 
equations we choose to make this adjustment in the NN "potential", rather than in the 
pion absorption/production channel. 

The necessity of eliminating double-counting from the theory is therefore forcing either: 

1. The modification of the amplitudes appearing within the equations, or 

2. The inclusion of explicit subtraction terms to compensate for the over-counting. 

If option 1 is pursued the result is that different amplitudes appear at different places 
within the theory: for instance, we have four nN amplitudes in the equation (|117|) . If option 
2 is taken then the presence of extra subtraction terms means that the driving term of the 
equation takes on a far more complicated form, since it now involves diagrams with more 
than one explicit pion. In either case the original goal of deriving coupled equations for the 
NN — 7rNN system which contain only one explicit pion and standard 1PI iViV and ttN 
amplitudes has been defeated. If such equations were to be derived, as was done by AM 
and AB, certain diagrams would need to be (incorrectly) included in two (or more) places in 
the theory. That these diagrams appear once and only once in the diagrammatic expansion 
is a direct consequence of the indeterminacy of the time-ordering of interaction vertices on 
different particles in a time-dependent field theory. 



VII. DERIVATION OF COUPLED EQUATIONS FOR THE NN - irNN SYSTEM 



Our next task is to obtain double-counting-free coupled integral equations for the am- 
plitudes which describe the processes: 



Ax + A 2 
(N 2 n) + Ax 
(A l7 r) + N 2 

(NtN 2 ) + 7T 



N± + N 2 
(N 2 n) + Ni 
(A l7 r) + N 2 

(NtN 2 ) + 7T 



(120) 



This will be achieved in two stages: firstly, integral equations governing the amplitudes 
which describe the reactions NN -> NN, NN -> NNir, NNn -> NN, NNn -> NNn will 
be found. Once this is accomplished the residues of these amplitudes at the appropriate 
poles will be taken, in order to derive equations connecting the two-fragment amplitudes for 
the above processes. 
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A. Deriving integral equations for the two and three-body amplitudes 

In this section the first of these two steps is performed: the derivation of integral equations 
for the two and three-body amplitudes. We have, of course, already begun this task by 
deriving equations for the amplitudes T^ N , T^ N and F^ . But, equations are now needed, 
not only for these amplitudes, but also for the connected NNn — > NNn amplitude 
the connected NN — > NNn amplitude F^ , and the connected NNn — > NN amplitude, 
FW. (Note that since all the amplitudes under consideration involve a two-nucleon initial 
state, by conservation of nucleon number they are automatically one-particle irreducible.) 



1. FN* and 



Firstly, consider the connected 1PI NN — > NNn amplitude F^ . Since the parameters 
of the Taylor method for this case are n = 3, m = 2 and r = 2 no double- counting can arise 
and so the classification-of-diagrams technique may be applied to this amplitude without 
adaptation. This process produces the equation: 



:i2D 



where the first term is generated by C±, the second by C*2, and the third by C5. Classes C3 
and C 4 are empty, since, as has already been observed above, all one-to-one amplitudes are 
zero. 

Substituting the expression obtained for F^ given in Eq. ( J3TD then produces: 
F (i)t = £ f V\j) dld2 T$l + \Y,t^\a)dl%dJ ai f^\i) 



a/3i 



[1 + d\d2T 1 ^ N ) 



(122) 



Taking the adjoint of Eq. ( |122| ) then gives an expression for the 1PI NNn — > NN amplitude, 
F^ l \ in terms of two-body t-matrices and dressed nNN vertices. 



2. MW 



In this case m = n = 3 and r = 2, and therefore Taylor's original method may again be 
applied, giving: 



M (l) = M (2) + F (2)t di4F (l) + F (^ dl d 2 



i=l 



+ 



d 1 d 2 FW, (123) 



where the four terms come from, respectively, C\, C2, C3 and C5. Note that C4 is empty. 
Note also that we have made use of Eq. ( |B6| ) for here. 

In order to derive an expression for AfW in terms of two-body t-matrices and dressed 
nNN vertices an equation for must be derived. In the case of the parameters 
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of the Taylor method are n = m = r = 3 and so one might expect that, since n = r, 
double-counting will occur in C4 and C 5 . However, since all of the particles are fully dressed 
no three-cut may intersect more than one line from the final state and so in every diagram 
in C4 and C5 we have Sf = 1. Consequently, the condition, n > 2s f is always satisfied and 
so it is guaranteed that no double-counting occurs even if Taylor's original method is used. 
When the original Taylor method is applied to the amplitude it yields: 



M (2) = M (3) + 



d\d 2 d 7T 



+ E* (1) (/?)V 



(c) 



(124) 



where the superscript ^ indicates that the amplitude resulting from the multiplication must 
be connected. 

Since we have consistently ignored three-body NNn forces throughout this calculation, 
we set to zero, thus producing the equation: 



M (2) =J2t ( - 2 \a)d- 1 d 1 d 2 dJ a pt {1 \f3)dp 1 + J2t {2) (a)d~ 1 d 1 d 2 d 7r M (2) . 

a/3 a 

This equation may then be formally solved, in order to yield: 

M (2) =Y j t^\a)d~ 1 d l d 2 d^d l d 2 dJ l \(3)d- p \ 

a/3 



(125) 



(126) 



where U^l are the covariant AGS amplitudes, which obey Eqs. ([62]) and fl6"3"|). 

Using Eqs. (|B"TD, (|122| ) and ( |126|) in Eq. (|123| ) an expression for M^ l > in terms of the irN 
and NN amplitudes and the dressed nNN vertex may be derived. 



B. Coupled equations for the NN — > NN and NNtt — > NN reactions 

Eqs. (0), ( |117| ), ( p.22| ) and ( |123| ) are equations for the two and three-body amplitudes of 
interest. However, our goal is to derive equations connecting the amplitudes for the collision 
of two fragments, rather than equations containing the three-body amplitudes. So consider 
firstly the reactions in which a two-fragment final state that includes a pion is produced 
from an initial iViV state: 

f (iV 2 7r) + Ni 

Ni + N 2 — > < (JVitt) + N 2 (127) 
{ (JViJV 2 ) + 7T 

Suppose that the final state of interest is one in which the particle A is a spectator while the 
other two particles are in a bound or resonance state. If the energy of this final state is E\ 
then the matrix element for such a process, X® N , is given by fl32| : 

*xn = {^)(<PW\(xx\dx 1 did 2 d 7T Res x ^d^d^ 1 F% Y \^ D ) (128) 

where [0(A)) is the form factor for the formation of the bound or resonance state; \x\) and 
d\ are, respectively, the wave function and dressed propagator for a free A particle and \ipn) 
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is a wave function for the initial pair N% and N 2 . Note that so far we have treated all 
particles as distinguishable, hence the initial state NN wave function must be the one for 
distinguishable particles, and the form-factor 1 0(A)) should not be anti-symmetrized, even if 
A represents the channel (NiN 2 ) + 7r. Furthermore, we indicate the fact that the equations 

for assume distinguishable particles by placing the subscript d on 

In order to derive an equation for X® N the residue of Eq. ( |122| ) must be taken at the 
final state A pole. It is therefore necessary to make explicit the analytic structure "near" 
this pole of each of the amplitudes which appear last in the various terms on the right-hand 
side of Eq. (|122|). To facilitate this we first write: 



if a = 1,2 



mt \ - J t(1) ( a ) - t {l \a)d N d^v x (a) - v x {a) i 
J {a) -\tW(a) if a = 3 



(129) 



in Eq. (|122|) . Now, if p a is the total four- momentum of the two-body sub-system labeled by 



a then it may be shown that the analytic structure of t^(a) "near" the a pole, p 2 a = M ( 



2 



is given by: 



tM{a) vl ^ Mi |0(a)) 2 % _ m2 (0(qO1 + (regular terms). (130) 

Pa a 



When the three-body operator ^^(a)^ 1 is considered it is found that if the four-momentum 
of the spectator a is k a , the a pole is at the point (P — k a ) 2 = M 2 where P is the total 
four- momentum of the system. Near that point: 



t (1 \a)d a \k a ) {P ka l~* M <* d a \k a )\0(a))\xa(k a )) ^ p _ - % _ M2 d a (k a )((l)(a)\(x a (k a )\d a \k a ) 
+ (regular terms). (131) 

The amplitudes other than which appear last in some term in Eq. ( ggg ) are v and 
. Since the former is the crossed term of the n — N potential and the latter is the dressed 
N — > ixN form factor, neither v x or have a pole at any of the points (P — k\) 2 = E\. 
(The exception to this statement occurs if t^ N includes a resonance with the same quantum 
numbers as the nucleon. In this case care must be exercised as the residue may receive a 
contribution from = (l + t^yG^d,,-)/^.) Consequently, any term in which a v x or 
is the last amplitude disappears upon the taking of the residue. It follows that Eqs. (|128| ) 
and (|122|) lead to the following form for X® N : 

*XN = ((f>W\(xx\d x 1 d 1 d 2 d n T° N \ij D ), (132) 

where the amplitude T® N is given by: 

T? N = £(! - v x {X)d^ 1 d 1 d 2 d 7r )8 Xj f^\j)d- l (l + d x d 2 T° N ) 
j 

+ E U^d 1 d 2 dJ^(P)d^% l did 7T f w \i)(l + dxd 2 T§ N ). (133) 

pi 

Here we have replaced the A^A^ scattering amplitude T^ N by the operator T$ N , since, as 
discussed above, the 1PI A^A^ scattering amplitude is the full A^A^ scattering amplitude. The 
AGS equations may now be used in order to obtain: 
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rpD 

XN 



X 



- v x (\)d~ x 1 d 1 d 2 d n )S X3 f^\j)dj 1 -Y.^v x {a)d- l d 3 dJ a] f^\ 3 ) 

3 &j 

1 + d 1 d 2 T^ N ) +Y,hat {1) {<x)da 1 didAT° N . 



(134) 



Note that this equation is equivalent to that derived by Avishai and Mizutani (AM) |TIJ 
and Afnan and Blankleider (AB) ||11|| , except for the presence of the term: 



Y / v x {a)d- 1 d 3 dJ aj f^\j) 



[1 + did 2 T® N ), 



(135) 



which has been included in order to remove the double-counting present when the equations 
derived by AM and AB are used in a time-dependent perturbation theory. 

The equation derived for T® N contains the amplitude T$ N . Therefore it is now necessary 
to derive an equation for T$ N . Recall that Eqs. ( |ITD and ( |117| ) are the results given by the 
modified Taylor method for T NN . Substituting ( |117| ) into ( |Tl~D yields: 



Tnn = (Vqpe -D 1 -D 2 -D n -X-B)(l + d^T^) + ( £ '(j')^^ '(a)C 1 

. jai 



did 2 d n 



1 - v x {a)d- l )d 1 d 2 dj ai f^\i)di l + U^d x d 2 dJ^\^%dJ pi f^\i) 

& 

(l + dxd 2 T^ N ). (136) 
Using Eq. (|133|) for T aN then reveals that: 

Tg N = V(l + d x d 2 T° N ) +J2f [1) U)S j ad- j dJ^(a)d- 1 d 1 d 2 d^ N , (137) 



J a 



where V is given by Eq. ( |119| ). 

Recall that in the discussion after Eq. ( |117| ) it was emphasized that the final result for V 
is dependent upon the treatment of the iViV «-> NNn sector of the theory. Similarly, since 
Eq. (|137|) results from a consistent treatment of the coupled iViV — nNN system one cannot 
simple-mindedly modify this result by, for instance, merely dropping the second term on 
the right-hand side. Any such change must be developed in a coherent way throughout the 
theory, rather than imposed in an ad hoc way on the final equations. 

Because we have used the modified Taylor method and so eliminated all double-counting 
the set of four equations given in Eqs. ( |134| ) and ( |137| ) are double-counting-free coupled 
four-dimensional integral equations for the iViV — ► tiNN and iViV — > iViV channels. 



C. Coupled equations for the ttNN — > ttNN and NN — ► NNir reactions 

We now turn our attention to the nNN —> uNN and iViV — > nNN channels, and 
attempt to derive coupled equations for the reactions: 
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(N 2 n) + Ni 
(N l7 r) + N 2 

(NlN 2 ) + 7T 



and: 



N t + N 2 



(N 2 n) + ^ 

(iViTT) + N 2 
(NiN 2 ) + 7T 



(iV 2 7T) + Nl 

(JVitt) + iV 2 

(JV X iV 2 ) + 7T. 



138) 



(139) 



The two-to-two amplitude for any of the reactions ( 138|) is found by taking the right and 
left residue of at the relevant poles. That is, the matrix element for a transition from 



a state in which particle A is a spectator to a state in which particle \i is a spectator is XjP,, 



given by: 



X?„ = (-i) 2 ((f)(X)\(xx\d x 1 did 2 d 7T Res x ^ i e d\d2d^d x l ,d lx 1 did 2 d 7T 



Ref V pole^^d/l^))!^)- 



(140) 



In order to derive an equation for TR, Eq. (|123|) is used to glean an expression for 
in terms of the input to the equations: t^N, ^jvjv and J^. Using the facts about the analytic 
structure of these amplitudes which were given above and taking right and left residues of 
this expression as per Eq. (|140|) , it is found that: 



K = {m^xx^did^T^d^d- 1 ^^), 



(141) 



with: 



rpD _ rr(2) 
J A/i — ^A/t 



2(1 - v x (X)d^d 1 d 2 d 7r )S Xj f (1)1 {j)dj 1 

d\d 2 



+ 2 uSd 1 d 2 dJ 1 \5)dJ 1 d- j dJ 5j f^\j) 



(T^dKh + 1) 



Y l f (1 \i)S i7 d I dJ^(j)d- 1 d 1 d 2 d^ 
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+ 2/ (1) W^ rl ^(l - dxd^i^d- 1 ) 



(142) 



We may also define Xj^j^ as the two-fragment matrix element for transition from the two- 
body state in which particle is a spectator to the state containing two (distinguishable) 
nucleons: 



(143) 



X Nn = {-i){^D\F { D } d^dxd^Res^ ^^d^d^^^x^d^ . 
Once more it is found that: 

X^ = (MT^d^d- 1 ^))^), (144) 
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this time with: 



T Nu - {T^ N did 2 + 1) 



Y,f {l \^^d 3 dJ^\a)d- a 1 d 1 d 2 d^ 



jot 



+ J2f {1) U)dJ 1 S j ^-d 1 dAv x ( f M)d; 1 ) 



(145) 



Using the covariant AGS equations and Eq. ( |145[ ), Eq. ( |142| ) may now be simplified to: 



-iD 
Njj, 



K = h^d^d; 1 + ]T(1 - v x (X)d^d 1 d 2 d^8 Xj f {1) \j)d- 1 d 1 d 2 Ti 

j 

-Y.^ X ^)d~ 1 d 3 dJ aj f^\ 3 )d 1 d 2 T^ l + J2Sx a t {1 \a)d- 1 d 1 d 2 d 7T T^. (146) 



a-3 



Note that once again, this equation is equivalent to that derived by AM [|H]] and AB 
ljj , but with terms subtracted in order to eliminate the double- counting. 

Finally, an equation for is needed in order to close the set of equations for the 
NN — 7rNN system. Beginning with Eq. ( |145| ) for T$ , we use the covariant AGS equations 
and Eq. (|137|) for T$ N in order to obtain: 



= E/ (1) 0')^>(1 - dtd^i^d' 1 ) + Vd x d 2 T L 
j 

+ E f {1 \j)^ a d 3 dj^\a)d- 1 d l d 2 d^ 



D 



(147) 



j ° 



with V given by Eq. ( |119| ). Note the changes in this equation as compared with AM's (HI 
and AB's IffTI result: 



1. ThetW (a) appearing before f^(j) has been replaced by a t^(a). 
2- Vope has been replaced by V. 

3. The term J2j f^ 1 Kj)^jfj,djd n v x (fi)d~ 1 has been subtracted from the driving term. 



VIII. ANTI-SYMMETRIZATION OF THE COUPLED EQUATIONS 

At this stage the amplitudes T$ N , T^ N , T$ and TR are unphysical, since they have 
not been anti-symmetrized. In this section we derive equations for the anti-symmetrized 
amplitudes. These anti-symmetrized equations apply to identical nucleons. 

In order to anti-symmetrize the equations we must first derive the correct anti- 
symmetrization procedure. It is by no means self-evident that the ordinary non-relativistic 
procedure of taking appropriate linear combinations of unsymmetrized amplitudes will still 
be valid for these amplitudes which are fully relativistic and obey four-dimensional inte- 
gral equations. The only way to discover the correct anti-symmetrization procedure is to 
examine the Feynman Rules for the amplitudes. Examination of the Feynman Rules for a 
two-fermion Green's function shows that if 
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G%M,p 2 ; P i,P2) (148) 

is the Green's function which is obtained when two nucleons are treated as distinguishable 
particles then: 

Gnn(p' 1 ,p' 2 ;PuP2) = G NN (p' 1 ,p' 2 ;p 1 ,p 2 ) - G% N (p' 2 ,p[;p 1 ,p 2 ). (149) 

Applying LSZ reduction to each side of this equation, in order to obtain an equation 
connecting the amplitudes involved, then gives: 

T NN = (1-P 12 )T° N , (150) 

where T$ N is the amplitude for distinguishable particles and P\ 2 acting on any state inter- 
changes the roles of particles one and two in that state. Using the fact that: 



;i-p 12 ) 2 = 2(i-p 12 ), (i5i) 



and that: 



[P 12 ,T^] = 0, (152) 



Eq. ( |150| ) gives: 



T NN = ~(1 - P 12 )Tg N (l - P X2 ). (153) 



A similar argument for the NN — > NNtt and NNw — > NNn Green's functions shows: 

F (1)t = ^(1-^2)^(1-^2) (154) 

M (1) = ^(1 - P 12 )M%\l - P 12 ). (155) 
Consequently, if we now write: 

^physical _ (0(2)1(^1^^8^^ ^d^F^ Vd>, (156) 



we find that: 



where: 



^physical = m) \ {x2 \ dA T AN ^)^ (157) 



Tan - -^{T 2 D N - Pi 2 T® N ) , (158) 



and is a fully anti-symmetrized NN wave-function: 



-^(I-PuMd)- (159) 
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(Note that the choice of final-state spectator in Eq. ( |156| ) makes no difference to the result, 
beyond an overall minus sign.) Similarly: 

^physical _ ^( 2 )|(x 2 Mi4Res (JVi7r) ^^d^M^d^Res^^ pole rfic/ 2 |0(3)D) \Xn) (160) 

= {^{{xMid^TAddtd^id))^), (161) 

where: 

T Ad = ^{Tg-P X2 T°), (162) 

= -L(l_p 12 )|^(3) D >. (163) 

We now take the definitions of the distinguishable particle matrix elements X 22 , X^ 2 , 
X 32 , X NN , X 33 , X 3N and X N3 , replacing everywhere the distinguishable particle ampli- 
tudes with their identical-particle counterparts ( |153| )-(| l~5l)| ). This yields definitions for the 

. j . • i , • i . , . ^physical ^physical ^-physical ^physical ^physical 

identical-particle matrix elements X AA J , X^ A J , X*[ A J , X NN J , J , 

X^^ S ^ ca ^ and x|)^ S * ca ^. When taken together with the definitions X^^ S ^ ca ^ and 

X]y^ S * Ca ^ given above this produces nine definitions which may be written in the form: 



where A,B = N, A, d, 



^physical = ^ a \ Tab \^ b) ^ (164) 



\*I>n) = Wnn) (165) 
|^a> = ^410(2)) |X2> (166) 
\iP d ) = d 1 d 2 \<P(d))\ X .); (167) 



and: 



rp rpD . rp rpD . rp rpD . rp rpD . / -| /?Q\ 

J-NN — I NNi l-Nd — J-N?,! J-dN — J- 3 Ni I dd — J- 33 , (.IDS J 

Tan = —^(T^ - P 12 T 1N ); T NA = —f=(T m - T N1 P 12 ); (169) 

T Arf = -J=(P£ - P 12 P^); T dA = ±={Tg - TiP 12 ); (170) 

Taa = Pg - TgP 12 . (171) 

Equations for these amplitudes may then be derived from the Eqs. ( |134| ), ( |137| ), (|145[) 
and (|146|) for P^, T$ N , T® and TP. We find that the set of coupled equations: 

T AN = V2 [f(l)d^ - Y] (1 + d x d 2 T NN ) - P 12 t{2)d 1 d„T AN + y/2t(Z)d x d 2 T dN , (172) 



P 



1 + d x d 2 T NN ) + v / 2t(2) ( ii^rA J v, (173) 



Pvjv = y(l + didaPjviv) + V2f(l)d 2 dJ^2)d 1 d n T AN + ]T /(k^d^d^T^, (174) 
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describe the coupled A A — > nNN and iViV — > A A two-fragment processes. Meanwhile the 
equations: 



AA 



-P 12 t(2)d 1 ^T AA + V2t(3)d 1 d 2 T dA 



d\d 2 Tj^ A 



E/ t (*K 1 -y 



did 2 T NA + V2t(2)did n T AA 



T NA = V2f(l)di\l - d 1 d w v x (2)) + Vd 1 d 2 T NA + y^2f{l)d 2 dJ^{2)d 1 d n T AA 

+ H f(k)dkdnt(3)did 2 T dA 
k 

T Ad = V^d^d^d- 1 + [^(l)^ 1 - ^] did 2 T M - Pi2t{2)d 1 d n T Ad 
+V2t{3)d 1 d 2 T dd 



T dd 
TNd 



k 



did 2 T Nd + V2t{2)did n T Ad 
J2 f(k)d^ + Vd 1 d 2 T m + V2f(l)d 2 dJ\2)d 1 d 7T T Ad 

k 

+ H f{k)dkd n t(3)did 2 T dd , 



describe the NNn — > NNn and A A — > NNn two-fragment processes. Here: 

y = 5>*(fcK/t(*). 



(175) 
(176) 

(177) 

(178) 
(179) 

(180) 
(181) 



Note that the t-matrices and vertices used as input to these equations are always the one- 
particle irreducible amplitudes, /W and . An apparent exception is the inclusion 
of the modified one-particle irreducible t-matrix ft 1 " , which is used in some places in the 
equations. But, is, of course, defined in terms of t^ 1 ' and /W by Eqs. (|64f) and (|59|) . 



IX. SPECIFICATION OF INPUT AMPLITUDES AND IMPLICIT INCLUSION 

OF OTHER HADRONIC STATES 



Equations (|172|) - (|180|) are a set of two-fragment scattering equations for the processes: 

A + A + ttI fN + N + iv , . 

N + N j— * \ N + N ( 182 > 

Each Feynamn diagram in which no more than one explicit pion appears is included once 
and only once in the equations. Hence these four-dimensional AA — nNN equations obey 
two and three-body unitarity and contain no double-counting, thus providing a framework 
in which predictions for scattering observables in the A A — nNN system can be made. 

The first step in making such predictions is the specification of the input amplitudes 
and In Appendix [B| the one-explicit-pion content of t^jj and is elucidated by the 
Taylor method, demonstrating that the ampitudes which need to be specified are, in fact, 
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and f^. Now, t^ N and could have their iu>o-explicit-pion content analyzed (see e.g., 
for the tcN case in time-ordered perturbation theory). However, even were we able to 



solve the resulting set of coupled ttN — tctcN equations together with Equations ( |172|) - (|180| ) 
the results would still only include the effects of pion exchange. It is well known that other 
effects, e.g. heavy-meson exchange and delta degrees of freedom, are important in nuclear 
forces, and it is not clear how such mechanisms would be included in this description. 

Therefore in order to implicitly include some of this additional physics we advocate 
calculating and t^ N in some other (hadronic or quark) model and/or parametrizing 
them in such a way as to simplify the solution of Eqs. (|172|) -( ISD ). Note that since, by 
definition, these amplitudes have no one-explicit-pion content, such an approach does not 
affect NN or irNN unitarity. Nor, if it is implemented correctly, should it introduce any 
double-counting. 

Firstly, the bare ttNN form factor may be taken from some model of QCD, as long 
as care is taken with any pionic corrections so that overcounting is avoided. Alternatively, 
may simply be parametrized. 

Secondly, two approaches may be taken with t^ N . On the one hand, a separable potential 

(2) 

may be postulated for and the free parameters of the separable form fitted to the ttN 
scattering data. (Note that it may be necessary to use irN scattering data in the P\\ 
channel in order to fix any parameters in f^.) In this way known mechanisms of the 
7r — N interaction, such as p and cu exchange and the A resonance, are implicitly included 
insofar as they contribute to the actual data. On the other hand, if we wish to know 
exactly what physics is included in the input tyN interaction a meson exchange model of 
txN scattering, with the A resonance built in, may be constructed (see e.g. |34|) and a 



separable expansion of the resulting amplitude made, via the Ernst-Shakin-Thaler technique 
|35|,36]. Such an expansion has already been accomplished for the three-dimensional Paris 



iViV potential |57 |. In the ttN case Pearce and Afnan have used the extension of Ernst, 



Shakin and Thaler's technique developed by Pearce |38[ in order to obtain a separable 
expansion of the tiN interaction in the Cloudy Bag Model [|59|] . The work of Rupp and 
Tjon shows that making a separable expansion of a Bethe-Salpeter amplitude is not 
computationally intensive, thus once a meson-exchange model of the nN system has been 
used to calculate scattering via the Bethe-Salpeter Equation a separable expansion of the 
resulting amplitude should not be a difficult task. 

Thirdly, so far we have regarded T^ N as part of the solution to Eqs. ( |172| )-( |1"8"0|) . How- 
ever, this amplitude is also needed in order to construct the kernel of these equations. Indeed, 
it is TxlfS appearance in the kernel (as t$ N ) which ensures that higher pionic processes, such 
as crossed two-pion exchange, are included in the solution to the coupled equations (see Fig- 
ure [Tjp. In theory it is possible to "bootstrap" this theory up and so generate an exact 
solution to the non-linear equation for T^ N . However even if this could be achieved, impor- 
tant physics would be missing from the resulting iViV interaction. In order to simultaneously 
include some of these missing mechanisms and snap the "bootstrap" , we advocate the same 
approach to the input t$ N as to 4-at- Either Eq. ([11]) for Tjp N should be used and a separable 
potential for T^ N fitted to the experimental data, or a four-dimensional covariant separable 
expansion of some meson-exchange NN amplitude should be made. As in the tcN case, 
since the input t^ N occurs only in the presence of a spectator pion such an approach will 
not destroy NN or NNir unitarity, nor will it, if carefully implemented, introduce double- 
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counting. The consistency of such a procedure can be checked by comparing the solution 
Tft N with the input %]y. 

If the amplitudes %]y> t%N an d are constructed in this way then (provided a separable 
expansion for v x is also made) Eqs. ( |1T2| ) ( pT80|) become a set of coupled Bethe-Salpeter 
equations for the iViV — tcNN system. This is completely analagous to the situation in the 



non-relativistic three-body Faddeev equations [IT|. By using separable expansions for the 
input amplitudes these may be reduced to the Lovelace equations IpEfl , which are a set of 
coupled two-body equations. 

One question which must be answered is which energy domain any separable expansion 
of t$ N and i£jy should be made in. Suppose that Figure [19] is a piece of a Feynman Diagram 



summed in the four-dimensional NN — nNN equations, with a the spectator particle label 
(a = n, N) and k a its four- momentum. The analytic expression for this Feynman diagram 
would therefore contain a piece: 

J d 4 k a ■ ■ ■ G & {P - k a ) t£\s a ) d Q (k a ) G & {P -k a )'-- (183) 

where 

P=W~8,$) (184) 

s a = (P-k a ) 2 = (V~s-k° a ) 2 -kl (185) 

is the total four-momentum in the three-body centre of mass, and Ga{P — k a ) is the product 
of the two propagators of the interacting particles. 

Since k^ E (— oo, oo) the two-body energy which is the argument of ranges from — oo 
to oo. However, suppose that a separable expansion for t^> has been made, i.e.: 

t\s a )=Y.9\^- 2 ^ (186) 

with gi an appropriate form factor implicitly dependent on the relative four-momentum of 
the two particles in the interacting subsystem. In order to perform the k Q integration we 
examine the analytic structure of the integrand in the complex plane. 

Apart from the dependence listed in Eq. ( |183| ), the only place k a appears in the integrand 



is via its presence in the form factors which occur immediately preceding and following the 
piece of diagram which is explicitly written in (|183|) . Unless these are nNN form factors 
their analytic structure in the relative momentum variable must give no contribution to 
unitarity, and therefore their /c^-dependence may be ignored. If one (or both) form factors 
is a ttNN form factor then there is a ttN cut (or cuts) to be added to the analytic structure 
in the lower half of the plane that is listed below. However, the existence of such cuts 
makes no difference to the final conclusion obtained here. Consequently, for the purposes of 
this argument, we ignore any fc Q -dependence of the form factors. 
The poles of the integrand in the lower half-plane are at: 



m 2 a + \k a \ 2 -ie (187) 



v/i+V^i + \k a \ 2 -ie (188) 
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while the two-body propagator contributes a cut beginning at: 



s + m a — ie (189) 

where m & is the mass of the two-body system. Provided that s is in the region for particle- 
particle scattering the piece of analytic structure giving the most significant contribution to 
the integral is the pole (|187|) . 



Thus the important two-body energies s a are those satisfying y/s^ < y/s - 
ment with the case of non-relativistic few-body physics |fL7| . Therefore if y/s is restricted 
to: 

2m < y/s < 2m + 2/i (190) 

then 7iN and iViV amplitude parametrizations based on information in the region below the 
respective pion production thresholds should give accurate results when used in the solution 
of Eqs. (|i7^)- ([T80|) . 

Finally, it is clear that NN interactions mediated by meson exchanges which are totally 
non-pionic are not included in our equations. These exchanges are known to be important 
effects in the medium to short range part of the nucleon-nucleon interaction. They could be 
included by extending the arguments given here for the pion in order to expose one-explicit- 
meson states for all mesons of importance. However, at this stage such an approach would 
be unnecessarily thorough, since we do not wish to make predictions for the reactions: 

NN + X <-> NN (191) 
X + X + d, (192) 

where X is any meson known to be important in one-boson-exchange iViV potentials! There- 

(2) 

fore we can afford to merely add a phenomenological heavy-meson exchange piece to the jy 
given by Eq. (|117|) , as was done in the three-dimensional NN — ttNN case by Avishai and 
Mizutani [jn],[£J. Note that some care must be taken in making this addition. For instance, 
the piece of p exchange arising from the exchange of two uncorrelated pions is separately 
included in the NN — irNN equations, therefore if a p is phenomenologically added to our 
formalism it will not necessarily have the same parameters as in NN OBE potentials. 



X. CONCLUSION 

In this paper we have derived coupled four-dimensional covariant scattering equations 
for the iViV — nNN system, using the modified Taylor method of classification of diagrams 
|21| . These equations sum all the covariant perturbation theory diagrams which include one 



explicit pion once and only once. They therefore: 

1. Are four-dimensional integral equations. 

2. Are covariant, not only on-shell, but also have off-shell covariance in the manner 
dictated by the Feynman graphical expansion. 

3. Are completely free of the double-counting problems of some previous four- dimensional 
NN - TriViV equations. 
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4. Obey NN and NNtt unitarity, by explicit construction. 

The double-counting subtractions found to be necessary are a consequence of attempting 
to derive a set of coupled equations which describe, in a unified way, the reactions: 

N + N ^ N + N 
N + N^N + N + tt 
N + N + n -> N + N + 7T, (193) 

by summing all covariant perturbation theory Feynman diagrams for the processes (|193| ) 
that involve one explicit pion or no explicit pions. Were one content merely to take the 
Bethe-Salpeter equation for iViV scattering and choose a set of Feynman diagrams for the 
kernel Tjp N no double-counting would arise. However, such an approach does not satisfy 
three-body unitarity, since it excludes certain s-channel two-particle irreducible diagrams 
which contain one explicit pion. Only by deriving a set of coupled equations which sum all 
covariant perturbation theory diagrams containing zero or one explicit pion(s), as we have 
done here, will three-body unitarity be obeyed. 

The modified Taylor method used to make this derivation is completely general and so 
could also be used to obtain equations for other few-hadron processes. For example, the 
ttN — tttcN problem and pion photoproduction would both be amenable to a derivation by 
this method. 

The double-counting exposed and eliminated here occurred in four- dimensional iViV — 
ttNN equations. The implication of this double-counting for the "standard" three- 
dimensional NN — ttNN equations is not entirely clear. On the one hand, an obvious 
way to obtain three-dimensional equations from the ones given here is to take the two and 
three-particle Green's functions appearing in these equations and replace them by three- 
dimensional Green's functions, as suggested by, for example, Blankenbecler and Sugar |22 



It is apparent that the equations thus produced will be different to the standard NN — ttNN 
equations, since they will still include subtractions for double-counting. On the other hand, 
the standard NN — -kNN equations may also be derived by classifying diagrams in time- 
ordered perturbation theory, as was mentioned in the Introduction. If this approach is taken 
the resulting equations do not double-count any time-ordered perturbation theory diagrams. 
Hence, one may question whether the subtractions which would appear were one to apply a 
three-dimensional reduction to our NN — nNN equations are really necessary. 

Subject to this issue of the relation to a three-dimensional calculation, our results may 
have ramifications beyond few-hadron processes. For instance, the double-counting present 
in the unsubtracted four-dimensional NN — nNN equations may have implications for pion 
absorption and scattering on larger nuclei, as has been discussed by, for example Kowalski 
et al. Pop . Therefore the remedies for this double-counting discussed above may well prove 



applicable in these larger systems. 

The set of equations obtained here are (apart from the minor point discussed in Sec- 
tion |VI B 1|) equivalent to those previously obtained by Kvinikhidze and Blankleider (KB) 



24] , |25[] , although we have anchored our derivation more firmly in Taylor's original work. 
Furthermore, the final set of coupled equations given here are ready for computation, being 
written in the form of two-fragment scattering equations. 

These scattering equations require as their input only the ir — N t-matrix and dressed 
nNN vertex. However, since the NN t-matrix is present in the equations both within the 
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kernel and as part of the solution it appears that some way of "bootstrapping" the theory up 
must be found. In Section |X| we discussed how this bootstrap might be avoided and came 
to the inclusion that making a separable expansion of a model NN amplitude at energies 
up to the pion-production threshold should provide a good input NN amplitude. Even once 
the bootstrap problem is resolved in this way the computation is still somewhat formidable, 
involving, as it does, the numerical solution of a set of coupled four-dimensional integral 
equations. Nevertheless, the equations derived here (together with those obtained by KB) 
represent the first complete and correct summation of the one-explicit-pion sector of the 
theory of the NN — irNN system. Not until these equations are solved will it be clear 
how much of the physics of the NN — tcNN system is attributable to the one-explicit-pion 
sector. Therefore, we believe that the numerical work involved here, although at first sight 
somewhat daunting, is necessary. 

Despite the fact that we have had to expose some NNirn states to remove double- 
counting on the NNn level we have, in general, refrained from entering the two-pion sector 
of the theory. This has meant that the amplitudes M^ 3 \ F^ 3 ^ and have been completely 
ignored. The modification to the above equations to include such "three-body" mechanisms 
is one obvious way in which these equations could be extended. 

We note that mesons other than the pion and baryonic states such as the A may be 
implicitly included in the equations via their (possibly parametrized) presence in input 
amplitudes. They may be included explicitly if the Lagrangian is extended to one involving 
all hadrons of interest. Such a modification of the underlying field theory merely results 
in a proliferation of possible diagrams, rather than any fundamental change in the way the 
derivation proceeds. 

However, such improvements are in the future. Until the equations as they stand here 
are solved, and it is seen how well they reproduce the experimental data in the NN — 
nNN system, we shall not know which, if any, of these more sophisticated mechanisms are 
necessary for a correct description of the NN — 7rNN system dynamics. The next stage of 
this work must therefore be an effort to solve the equations given here numerically. 

As a first step towards this goal we have recently used the formalism developed here to 
perform a calculation of <f>(p scattering in a scalar <p 2 a field theory. Numerical methods and 



results will be detailed in a later paper |43 



APPENDIX A: PROPAGATOR DRESSING 

In this appendix we discuss how to accomplish the dressing of all the particles involved. 
The procedure used is similar to that developed by Afnan and Blankleider [|Il"| , |T7|1 . 

Consider the amplitude A n ^ m for the m — > n transition. The diagrams contributing to 
A may be of two types: 

1. Those for which a self-energy contribution on one leg, which we consider to be leg i in 
the initial state, may be isolated. (The argument is exactly the same for self-energy 
contributions on legs in the final state.) These diagrams are of the general form: 

A n ^ m d\ Ej, (Al) 
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where d^ and Ej are the undressed propagator and self-energy of particle i. (If there 
are j nucleons in the initial state then for i < j d^ and Ej represent a nucleon prop- 
agator and self-energy, otherwise they represent a pion propagator and self-energy) 
See Figure ^0| for an example of such a diagram. 

2. The diagrams which cannot have a self-energy contribution isolated in this way. We 
denote these by: 

A n ^ m& (A2) 



Therefore, 



A n <— m — A n ^ m d^ Sj + A. n< — m fi\, (A3) 
=^ A n< — m d\ ^ = A n ^ m fydi. (A4) 



Since, 



d t = df\l - E^f)" 1 = (4 0) 1 - E,)- 1 = Zidf, (A5) 

where Zi is the wave function renormalization for particle i, and df is the renormalized 
propagator with unit residue at the pole of dj which corresponds to the physical mass of 
particle i. 

Repeating this procedure for all external legs leads to the following result: 
where: 



k 



G fe = n^ 5 , (A7) 

i=\ 

Gf = V,df\ (A8) 

and An^fn has had all the external bubbles removed from it, and factors of Z^ included in 
it, in order to make it agree with Eq. 

This result allows us to consider all initial and final state legs as fully-dressed, provided 
that we work with the amplitude A^^ which has no bubbles on these initial or final state 
legs. Bubbles may also be eliminated from the internal legs of An<-m by a similar procedure, 
as follows. Consider any set of diagrams contributing to A^-t-m which are (k — l)-particle 
irreducible but admit an internal fc-particle cut. According to the last internal cut lemma 
these diagrams may be expressed uniquely as: 

iWSfiV*. (A9) 

Using the same arguments as above any bubbles appearing in the initial state of A^k m ay 
be amputated and placed in G k . Consequently, we obtain for this sum: 

A^~ k G k A~ k ^. (A10) 

Therefore we may always work with fully-dressed particles, provided that we consider all 
amplitudes to be of the type An<-m- In this paper we have adopted this approach, and so 
all amplitudes are of the Ai^-m type, even though the tildes are never displayed explicitly. 
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APPENDIX B: INPUT TO THE NN - ttNN EQUATIONS: THE vriViV VERTEX 

AND THE vr - N AMPLITUDE. 



In this appendix we consider the two amplitudes which are required as input to the 
NN — nNN equations: the one-particle irreducible 7riViV vertex, and the one-particle 
irreducible tt — N t-matrix, t^. 

Firstly, consider Our aim is to apply the Taylor method to this vertex in order to 
derive an integral equation for it. Since we are examining the two-cut structure of we 
have m — 2,n — \ and r = 2. Consequently, since we are dealing with fully-dressed particles 
all one-to-one amplitudes are zero and so classes C 3 , C 4 and C 5 are all empty. Therefore 
double-counting in classes C4 and C5 cannot arise, even though n < r. The sum of class C\ 
in this case is clearly the two-particle irreducible ttNN vertex f( 2 \ The sum of class C2 is 
found, via the last internal cut lemma, to be: 

f {2) d N d^ ] N . (Bl) 

Therefore, Taylor's method gives the following equation for f^: 

/( D = /(2 ) + f m dNd ^l ( B2) 

Taylor's method could now be used in order to extract the structure of f( 2 \ but we prefer 
to not consider three-body unitarity for this amplitude. Instead the expression for will 
be extracted from the Lagrangian under consideration. 

However, we do wish to discuss the two-particle cut structure of the 1PI rt — N amplitude, 
t^, using the Taylor method. Observe that classes C3, C4 and C5 are again empty, since 
all particles are fully dressed. Applying the Taylor method, it is clear that the sum of class 

( 2) 

d is the 2PI amplitude iffr. Also, the LICL may be applied to class C 2 in order to obtain: 

t^jfdNd n t^lf; (B3) 
and therefore the following Bethe-Salpeter type equation for tfy is obtained: 

t*N = tnN + tnN d Nd n t^l- (B4) 

As mentioned above in the case of f( 2 \ the structure of t^ N may be investigated using similar 

techniques as have been used for t^ N . On the other hand, if we are not concerned with nnN 

(2) 

unitarity in the ir — N amplitude we may just extract some model-dependent result for 

from the Lagrangian under consideration. 

(i) 

ttTV 

Taylor's method to t^h in the same fashion as above, reveals that: 

fi = *S + / (1) W (0) - (B5) 

A brief examination of the zero-particle irreducible irNN amplitude, /'- -', shows that the 
only Taylor class of diagrams contributing to /(°) which is non-empty is C±, which sums to 
give the 1PI nNN amplitude, Therefore, 



The amplitude t^h is, however, only part of the full ir — N amplitude, t^. Applying 
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/(°) = /«. (B6) 
Consequently, the full tc — N amplitude ^Sr i s given by: 

« = S + / (1) W (1) , (B7) 

where tfy is the non-pole part of the it — N t-matrix, given by Eq. ([B3]), and d^f^ is 
the pole part of the n — N t-matrix. 



APPENDIX C: DETAILS OF THE CALCULATION OF a n C^ N2 ' } 

The portion of t^ N which is 2PR in the N' <— n' + n + N channel may be written, using 
the reverse of the LICL, the "first internal cut lemma", as: 

fgdid^pi (CI) 

where P^ is an iW 2 — > Nni^ amplitude which is two-particle irreducible in the s-channel 
and the channel: 

7V[ + N' <- N + 7T 2 + 7T 2 , 

as well as being 1PI in the channel: 

7T 2 + N «- N' + 7T 2 + 7l[. 

Consequently, we might expect the diagrams: 

f^{l)d^{2)d 1 d Vl P\2)d 2 d w J^{l) (C2) 

and: 

/«(l)^J«(2)d 1 ^ 1 Pt(2)d 1 ^ 2 ttV(2)rf 2 ^J« t (l) (C3) 

to be included in c\ fl C5 . Consider first (|C2|) . If this term is to be part of C\ the 
pion absorbed on N2' must be "hidden" in another amplitude, since if the diagram is to 
contribute to C5 it must take the form: 

f^(2)d 1 d n F^\ (C4) 
There are three ways of "hiding" the pion: 

1. Place it in an N — N t-matrix. Taking this approach shows that the diagram: 

/( 1 )(l)^tf ) (2)rf 1 tS(2)rf 7r ,rf 1 /W t (2)^d 2 /« t (l) (C5) 

is in Ci n CI , since it occurs in both c\ and the term: 

/W(2)444 1 ^i/ (1)t (2) (C6) 



ofC 5 



{N2 1 } 
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2. The second way of hiding the second pion is to place it in the it — N t-matrix. That 
is, to note that the diagrams: 

fg\l)d^d a f^2)d n d 1 ^d n fff\2)d n fS\l) and 
is included in both C\ and: 

/«(2)d 7r d 2 t«(l)4/ (1)t (2), (C8) 

which is part of C\ N2 \ 

3. Finally, one can "hide" the pion in a three-body force. This suggests that if the term: 

f^{2)d x d 1 Mi 3) d 2 dJ^\l) (C9) 
is included in C\ N2 ^ then a number of diagrams in c\ will also occur in C\ N2 ^ . Precisely 



which diagrams from c\ are included in (|C9f) though? The additional irreducibility 

(3) ~ 

constraints placed on M\ mean that only that part of c±, Ci, which is given by: 

~ Cl = /W(l)i 2 d 2 /«(2)d 1 i 1 pt(2)4 2 f t (l), (CIO) 

where P^ is 1PI in both the channels: 

N' + N <- n[ + n' 2 + tt 2 and tc[ + N ^ N' + n' 2 + tt 2 , 

is included in flC9Q . The contributions to C\ produced by the parts of which are 
1PR in these two channels must be dealt with separately. This, in fact, is the reason 
why possibilities 1 and 2 had to be dealt with above. Once this is done, however, only 
the part c\ of c\ is left for consideration. But, c\ C c\ H C\ if and only if ( |C9| ) is 
included in CI . Consequently, in this calculation, in which M{ is set to zero, the 
diagrams C\ are not in c\ fl C§ N2 \ 



Turning now to (|G3|) , similar arguments to the above show that the expression (|C3| ) is 



only part of c\ N2 ^ if a pure three-body force is included in the calculation. Therefore, for 
the purposes of this argument ( |C3|) is not in C$ N2 (1 c± . 

APPENDIX D: DETAILS OF THE CALCULATION OF c 2 n C 5 {7V2 } 

The first internal cut lemma may be used to show that the diagrams which contribute 
to fjvh and are 2PR in the NV <- Nl + N2 + A^2'-channel sum to give: 

f^ l \2)d 1 d^\ (Dl) 
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where is two-particle irreducible in the s-channel and the channel: 

tt' + ATI' <- Nl + N2 + N2', 
as well as being 1PI in the channel: 

N2' + N2^ Nl' + n' + Nl. 
The question now is, which portions of the two expressions: 

f^\l)d 2 d n jW(2)d 2 d n2 F^d 2 fif(l) (D2) 

and: 

/«(l)d 2 4J«(2)d 2 4 2 F( 2 ) t d 1 44 1 l r d 2 /W t (l) (D3) 

are also in C^ N2 ^? As in Appendix we argue that we are not interested in those portions 
of expressions ( |L)2| ) and ( |D3| ) which may be reexpressed in the form ( |C9| ) . This line of 
reasoning leads to the conclusion that the diagrams contained in ( |D2| ) and C§ N2 ^ are: 

4:H2)4ji:Hi)^C(i)rf 2 ^j(: )t (i)rf 1 ^ 1 /^ )t (2), (D4) 

and: 

f^{2)d l d w jS{l)d, 2 d 1 f^ (D5) 

where T^lf is not only one-particle irreducible in the s-channel, but also 1PI in the t-channel. 

Similarly, there are no diagrams which are in both ODD and Ci N2>} provided, once again, 
that no pure three-body force is included in the calculation. 
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FIGURES 



FIG. 1. The term on the right is included in the current NN — ttNN equations, while the term 
on the left, known as the Jennings term, is not included. 



FIG. 2. One covariant perturbation theory diagram which is double-counted if the original 

(2) 

Taylor method is used to construct an equation for with the two cuts which place it in both 

Cf 1} and Cf 2} . 



FIG. 3. If the two cuts shown are three-cuts then they will place certain diagrams which should 
only be included in C4 in C5 as well. 



50 



(2) 

FIG. 4. One possible source of double-counting: if i^y is one-particle reducible in the u-channel 
then the cut shown places this diagram in cj 71 ^ N2 \ Since the diagram arose as part of c\ N2 ^ N1 ^ 
such u-channel one-particle reducibility leads to double-counting. 



(2) 

FIG. 5. Another possible source of double-counting: if is one-particle reducible in the 
t-channel then the cut shown places this diagram in C^ N2 ^ Nl \ Since the diagram arose as part 
of sub-class c\ N1 ^ N2 ^ such t-channel one-particle reducibility would lead to double-counting. 



FIG. 6. An example of a diagram which would lead to double-counting of the type depicted in 
Figure || 
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FIG. 7. The first term in the sub-class c\ N1 ^, c\, with the cut which may lead to dou- 
ble-counting indicated. 



FIG. 8. The lowest order diagrams which contribute to ir — ir scattering in a Lagrangian with 
a irNN vertex and a ttN contact term. 



FIG. 9. The second term of c\ Nl \ C2, with the cut which may lead to double-counting indi- 
cated. 
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FIG. 10. The third term in the sub-class c\ Nl \ C3, with the cuts which may lead to dou- 
ble-counting indicated. 



FIG. 11. The fourth term in c\ N1 \ C4, and the cut which may lead to double-counting. 



FIG. 12. The two possible cuts which may lead to double-counting fifth term of C§ N1 ^ , C5. 
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(N2'\ 

FIG. 13. The two possible cuts which may lead to double-counting in the sixth term of C5 , 

c 6 . 



FIG. 14. The definition of the new vertex, /W* 
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FIG. 15. Two diagrams which are part of T^j^. They will both contribute to the dressing of 
the vertices in one-pion exchange when the one-pion exchange piece of T^ N is inserted, as shown 
here. Hence both diagrams will lead to double-counting, unless the vertex used in the one-pion 

(2) 

exchange piece of T^-Jv is adjusted. 
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FIG. 16. On the right we show the three terms D\, D2 and D^. On the left we depict, in each 

(2) 

case, the two different places they appear in the pion absorption/production piece of T^jy. 
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FIG. 17. In the lower half of this diagram we show the two terms B and X, while the upper 
half displays the two places they appear in the pion absorption/production piece of T^ N . 



FIG. 18. A diagram containing the input NN amplitude t$ N . Note that if the one-pion ex- 
change piece of the NN amplitude is considered this diagram yields the crossed two-pion exchange 
graph. 
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FIG. 19. A piece of an arbitrary Feynman diagram to be included in the four-dimensional 
NN - itNN equations. 



FIG. 20. An example of a diagram in which a self-energy contribution on one of the legs 
belonging to the initial state may be isolated. 
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